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ABSTRACT.  Since  its  launch  in  1999,  the  Far  Ultraviolet  Spectroscopic  Explorer  (FUSE)  has  made  over 
4900  observations  of  some  2500  individual  targets.  The  data  are  reduced  by  the  principal  investigator  team  at 
the  Johns  Hopkins  University  and  archived  at  the  Multimission  Archive  at  STScI  (MAST).  The  data  reduction 
software  package,  called  CalFUSE,  has  evolved  considerably  over  the  lifetime  of  the  mission.  The  entire  FUSE 
data  set  has  recently  been  reprocessed  with  CalFUSE  version  3.2,  the  latest  version  of  this  software.  This  paper 
describes  CalFUSE  version  3.2,  the  instrument  calibrations  on  which  it  is  based,  and  the  format  of  the  resulting 
calibrated  data  files. 


1.  INTRODUCTION 

The  Far  Ultraviolet  Spectroscopic  Explorer  (FUSE)  is  a 
high-resolution,  far-ultraviolet  spectrometer  operating  in  the 
905-1187  A  wavelength  range.  FUSE  was  launched  in  1999 
on  a  Delta  II  rocket  into  a  nearly  circular,  low-Earth  orbit  with 
an  inclination  of  25°  to  the  equator  and  an  approximately 
100  minute  orbital  period.  Data  obtained  with  the  instrument 
are  reduced  by  the  principal  investigator  team  at  the  Johns 
Hopkins  University  using  a  suite  of  computer  programs  called 
CalFUSE.  Both  raw  and  processed  data  files  are  deposited  in 
the  Multimission  Archive  at  STScI  (MAST).13 

CalFUSE  evolved  considerably  in  the  years  following 
launch,  as  our  increasing  knowledge  of  the  spectrograph’s  per¬ 
formance  allowed  us  to  correct  the  data  for  more  and  more 
instrumental  effects.  The  program  eventually  became  unwieldy, 
and  in  2002  we  began  a  project  to  rewrite  the  code,  incorpo¬ 
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rating  our  new  understanding  of  the  instrument  and  best  prac¬ 
tices  for  data  reduction.  The  result  is  CalFUSE  version  3,  which 
produces  a  higher  quality  of  calibrated  data  while  running 
10  times  faster  than  previous  versions.  The  entire  FUSE  archive 
has  recently  been  reprocessed  with  CalFUSE  version  3.2;  we 
expect  this  to  be  the  final  calibration  of  these  data. 

In  this  paper,  we  describe  CalFUSE  version  3.2.0  and  its 
calibrated  data  products.  Because  this  document  is  meant  to 
serve  as  a  resource  for  researchers  analyzing  archival  FUSE 
spectra,  we  emphasize  the  interpretation  of  processed  data  files 
obtained  from  MAST  rather  than  the  details  of  designing  or 
running  the  pipeline.  An  overview  of  the  FUSE  instrument  is 
provided  in  §  2,  and  an  overview  of  the  pipeline  in  §  3.  Sec¬ 
tion  4  presents  a  detailed  description  of  the  pipeline  modules 
and  their  subroutines.  The  FUSE  wavelength  and  flux  calibra¬ 
tion  are  discussed  in  §  5,  and  a  few  additional  topics  are  con¬ 
sidered  in  §  6.  A  detailed  description  of  the  various  file  formats 
employed  by  CalFUSE  is  presented  in  the  Appendix. 

Additional  documentation  available  from  MAST  includes 
the  CalFUSE  home  page,14  the  CalFUSE  Pipeline  Reference 
Guide,15  the  FUSE  Instrument  and  Data  Handbook,16  and  the 
FUSE  Data  Analysis  Cookbook.17 

2.  THE  FUSE  INSTRUMENT 

FUSE  consists  of  four  coaligned  prime-focus  telescopes, 
each  with  its  own  Rowland  spectrograph  (Fig.  1).  Two  of  the 
four  channels  employ  Al+LiF  optical  coatings  and  record  spec¬ 
tra  over  the  wavelength  range  ~990— 1 1 87  A,  while  the  other 


14  See  http://archive.stsci.edu/fuse/calfuse.html. 

15  See  http://archive.stsci.edu/fuse/pipeline.html. 

16  See  http://archive.stsci.edu/fuse/dhbook.html. 

17  See  http://archive.stsci.edu/fuse/cookbook.html. 
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Fig.  2. — FUSE  apertures  projected  onto  the  sky.  In  the  FPA  coordinate 
system,  the  LWRS,  HIRS,  and  MDRS  apertures  are  centered  at  Y  =  —  1 18.07", 
—  10.27",  and  +90.18",  respectively.  The  reference  point  (RFPT)  at  X  = 
+  55.18"  is  not  an  aperture;  when  a  target  is  placed  at  this  location,  the  three 
apertures  sample  the  background  sky.  With  north  on  top  and  east  on  the  left, 
this  diagram  corresponds  to  an  aperture  position  angle  of  0°.  Positive  aperture 
position  angles  correspond  to  a  counterclockwise  rotation  of  the  spacecraft 
about  the  target  aperture.  This  diagram  represents  only  a  portion  of  the  FPA; 
its  active  area  is  19'  x  19'. 


Fig.  1. — Schematic  of  the  FUSE  instrument  optical  system.  The  telescope 
focal  lengths  are  2245  mm,  and  the  Rowland  circle  diameters  are  1652  mm. 
(Figure  from  Moos  et  al.  2000.) 

two  use  SiC  coatings,  which  provide  reflectivity  to  wavelengths 
below  the  Lyman  limit.  The  four  channels  overlap  between  990 
and  1070  A.  Spectral  resolution  is  roughly  20,000  (X/AX)  for 
point  sources.  For  a  complete  description  of  FUSE ,  see  Moos 
et  al.  (2000)  and  Sahnow  et  al.  (2000b). 

At  the  prime  focus  of  each  mirror  lies  a  focal  plane  assembly 
(FPA;  Fig.  2)  containing  three  spectrograph  entrance  apertures: 
the  low-resolution  aperture  (LWRS;  30"  x  30"),  used  for  most 
observations,  the  medium-resolution  aperture  (MDRS;  4"  x  20"), 
and  the  high-resolution  aperture  (HIRS;  1.25"  x  20").  The  ref¬ 
erence  point  (RFPT)  is  not  an  aperture;  when  a  target  is  placed 
at  this  location,  the  three  apertures  sample  the  background  sky. 
For  a  particular  exposure,  the  FITS  file  header  keywords 
RA_TARG  and  DEC_TARG  contain  the  J2000.0  coordinates 
of  the  aperture  (or  RFPT)  listed  in  the  APERTURE  keyword, 
while  the  keyword  APER_PA  contains  the  position  angle  of 
the  —  Y  axis  (in  the  FPA  coordinate  system;  see  Fig.  2),  cor¬ 
responding  to  a  counterclockwise  rotation  of  the  spacecraft 
about  the  target  (and  thus  about  the  center  of  the  target 
aperture). 

The  spectra  from  the  four  instrument  channels  are  imaged 


onto  two  photon-counting  microchannel-plate  (MCP)  detectors, 
labeled  1  and  2,  with  a  LiF  spectrum  and  a  SiC  spectrum  on 
each  (Fig.  1).  Each  detector  is  comprised  of  two  MCP  segments, 
labeled  A  and  B.  Raw  science  data  from  each  detector  segment 
are  stored  in  a  separate  data  file;  an  exposure  thus  yields  four 
raw  data  files,  labeled  1A,  IB,  2A,  and  2B.  Because  the  three 
apertures  are  open  to  the  sky  at  all  times,  the  LiF  and  SiC 
channels  each  generate  three  spectra,  one  from  each  aperture. 
In  most  cases,  the  nontarget  apertures  are  empty  and  sample 
the  background  sky.  Figure  3  presents  a  fully  corrected  image 
of  detector  1A  obtained  during  a  bright-Earth  observation.  The 
emission  features  in  all  three  apertures  are  geocoronal.  Note 
that  the  LiF  1  wavelength  scale  increases  to  the  right,  while 
the  SiC  1  scale  increases  to  the  left.  The  Ly/3  X1026  airglow 
feature  is  prominent  in  each  aperture. 

Two  observing  modes  are  available:  in  photon-address 
mode,  also  known  as  time-tag  or  TTAG  mode,  the  X  and  Y 
coordinates  and  pulse  height  (§  4.3.7)  of  each  detected  photon 
are  stored  in  a  photon-event  list.  A  time  stamp  is  inserted  into 
the  data  stream,  typically  once  per  second.  Data  from  the 
entire  active  area  of  the  detector  are  recorded.  Observing 
bright  targets  in  time-tag  mode  can  rapidly  fill  the  spacecraft 
recorder.  Consequently,  when  a  target  is  expected  to  generate 
more  than  ~2500  counts  s  1  across  all  four  detector  segments, 
the  data  are  stored  in  spectral-image  mode,  also  called  histo- 
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Fig.  3. — Image  of  detector  segment  1A  during  a  bright-Earth  observation.  All  lines  are  geocoronal.  Note  the  strong  Ly/3  (1026  A)  feature  in  each  spectrum. 
The  data  have  been  fully  corrected  for  detector  and  other  distortions.  Extended-source  extraction  windows  for  all  three  apertures  in  both  the  LiF  and  SiC  channels 
are  marked;  point-source  extraction  windows  are  somewhat  narrower  in  Y.  Instrumental  astigmatism  is  responsible  for  the  bow-tie  shape  of  each  spectrum.  The 
region  shown  corresponds  to  detector  pixels  900-15,300  in  X  and  0-915  in  Y. 


gram  or  HIST  mode.  To  conserve  memory,  histogram  data  are 
(usually)  binned  by  eight  pixels  in  Y  (the  spatial  dimension), 
but  unbinned  in  X  (the  dispersion  dimension).  Only  data  ob¬ 
tained  through  the  target  aperture  are  recorded.  Individual  pho¬ 
ton  arrival  time  and  pulse  height  information  is  lost.  The  orbital 
velocity  of  the  FUSE  spacecraft  is  7.5  km  s  !.  Since  Doppler 
compensation  is  not  performed  by  the  detector  electronics,  his¬ 
togram  exposures  must  be  kept  short  to  preserve  spectral  res¬ 
olution;  a  typical  histogram  exposure  is  about  500  s  in  length. 

The  front  surfaces  of  the  FPAs  are  reflective  in  visible  light. 
On  the  two  LiF  channels,  light  not  passing  through  the  apertures 
is  reflected  into  a  visible-light  CCD  camera.  Images  of  stars 
in  the  field  of  view  around  the  apertures  are  used  for  acquisition 
and  guiding  by  this  camera  system,  called  the  Fine  Error  Sensor 
(FES).  FUSE  carries  two  redundant  FES  cameras,  which  were 
provided  by  the  Canadian  Space  Agency.  FES  A  views  the 
FPA  on  the  LiF  1  channel,  and  FES  B  views  the  LiF  2  FPA. 
During  initial  checkout,  FES  A  was  designated  the  default  cam¬ 
era  and  was  used  for  all  science  observations  until  it  began  to 
malfunction  in  2005.  In  July  of  that  year,  FES  B  was  made  the 
default  guide  camera.  Implications  of  the  switch  from  FES  A 
to  FES  B  are  discussed  in  §  6.1. 


3.  OVERVIEW  OF  CALFUSE 

The  new  CalFUSE  pipeline  was  designed  with  three  prin¬ 
ciples  mind:  the  first  was  that  to  the  extent  possible,  we  follow 
the  path  of  a  photon  backward  through  the  instrument,  cor¬ 
recting  for  the  instrumental  effects  introduced  in  each  step.  The 
principal  steps  in  this  path,  together  with  the  effects  imparted 
by  each,  are  listed  below.  Most  of  the  optical  and  electronic 
components  in  this  list  are  labeled  in  Figure  1. 

1.  Satellite  motion  imparts  a  Doppler  shift. 

2.  Satellite  pointing  instabilities  shift  the  target  image  within 
(or  out  of)  the  aperture. 

3.  Thermally  induced  mirror  motions  shift  the  target  image 
within  (or  out  of)  the  aperture. 

4.  FPA  offsets  shift  the  spectrum  on  the  detector. 

5.  Thermally  induced  motions  of  the  spectrograph  gratings 
shift  the  target  image  within  (or  out  of)  the  aperture. 

6.  Ion-repelling  wire  grids  can  cast  shadows,  called  “worms.” 

7.  Detector  effects  include  quantum  efficiency,  flat  field,  dead 
spots,  and  background. 

8.  The  spectra  are  distorted  by  temperature,  count-rate,  time, 
and  pulse-height  dependent  errors  in  the  photons’  measured  X 
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and  Y  coordinates,  as  well  as  smaller  scale  geometric  distortions 
in  the  detector  image. 

9.  Count-rate  limitations  in  the  detector  electronics  and  the 
Instrument  Data  System  (IDS)  data  bus  are  sources  of  dead 
time. 

To  correct  for  these  effects,  we  begin  at  the  bottom  of  the  list 
and  (to  the  extent  possible)  work  backward.  First,  we  adjust  the 
photon  weights  to  account  for  data  lost  to  dead  time  (step  9) 
and  correct  the  photons’  X  and  Y  coordinates  for  a  variety  of 
detector  distortions  (step  8).  Second,  we  identify  periods  of  un¬ 
reliable,  contaminated,  or  missing  data.  Third,  we  correct  the 
photons’  X  and  Y  coordinates  for  grating  (step  5),  FPA  (step  4), 
mirror  (step  3),  and  spacecraft  (step  2)  motions.  Fourth,  we 
assign  a  wavelength  to  each  photon  based  on  its  corrected  X 
and  Y  coordinates  (step  5),  then  convert  to  a  heliocentric  wave¬ 
length  scale  (step  1).  Finally,  we  correct  for  detector  dead  spots 
(step  7),  model  and  subtract  the  detector  and  scattered-light 
backgrounds  (step  7),  and  extract  (using  optimal  extraction,  if 
possible),  flux  calibrate  (step  7),  and  write  to  separate  FITS 
files  the  target’s  LiF  and  SiC  spectra.  Note  that  we  cannot 
correct  for  the  effects  of  worms  (step  6)  or  the  detector  flat 
field  (step  7). 

Our  second  principle  was  to  make  the  pipeline  as  modular 
as  possible.  CalFUSE  is  written  in  the  C  programming  language 
and  runs  on  Solaris,  Linux,  and  Mac  OS  X  (vers.  10.2  and 
higher)  operating  systems.  The  pipeline  consists  of  a  series  of 
modules  called  by  a  shell  script.  Individual  modules  can  be 
executed  from  the  command  line.  Each  performs  a  set  of  related 
corrections  (screen  data,  remove  motions,  etc.)  by  calling  a 
series  of  subroutines. 

Our  third  principle  was  to  maintain  the  data  as  a  photon  list 
(called  an  intermediate  data  file,  or  IDF)  until  the  final  module 
of  the  pipeline.  Input  arrays  are  read  from  the  IDF  at  the  be¬ 
ginning  of  each  module,  and  output  arrays  are  written  at  the 
end.  Bad  photons  are  flagged  but  not  discarded,  so  the  user 
can  examine,  filter,  and  combine  processed  data  files  without 
rerunning  the  pipeline.  Like  all  FUSE  data,  IDFs  are  stored  as 
FITS  files  (Hanisch  et  al.  2001).  The  various  file  formats  em¬ 
ployed  by  CalFUSE  are  described  in  the  Appendix. 

A  FUSE  observation  consists  of  a  set  of  exposures  obtained 
with  a  particular  target  in  a  particular  aperture  on  a  particular 
date.  Each  exposure  generates  four  raw  data  files,  one  per  de¬ 
tector  segment,  and  each  raw  data  file  yields  a  pair  of  calibrated 
spectra  (LiF  and  SiC),  for  a  total  of  eight  calibrated  spectral 
files  per  exposure.  Each  raw  data  file  is  processed  individually 
by  the  pipeline.  Error  and  status  messages  are  written  to  a  trailer 
file  (described  in  §  4.10).  Spectra  are  extracted  only  for  the 
target  aperture  and  are  binned  in  wavelength.  Binning  can  be 
set  by  the  user,  but  the  default  is  0.013  A,  which  corresponds 
to  about  2  detector  pixels,  or  one-fourth  of  a  point-source  res¬ 
olution  element.  After  processing,  additional  software  is  used 
to  generate  a  set  of  observation-level  spectral  files,  the  ALL, 
ANO,  and  NVO  files  described  in  §  4.11.  A  complete  list  of 


FUSE  data  files  and  file-naming  conventions  can  be  found  in 
the  FUSE  Instrument  and  Data  Handbook.  All  of  the  exposures 
that  constitute  an  observation  are  processed  and  archived 
together. 

Investigators  who  wish  to  reprocess  their  data  can  retrieve 
the  CalFUSE  source  code  and  all  associated  calibration  files 
from  the  CalFUSE  home  page.  Instructions  for  running  the 
pipeline,  and  detailed  descriptions  of  the  calibration  files,  are 
provided  in  the  CalFUSE  Pipeline  Reference  Guide.  Note  that 
within  the  CalFUSE  software  distribution,  all  of  the  calibration 
files,  including  the  FUSE.TLE  file  (§  4.2),  are  stored  in  the 
directory  v3.2/calfiles,  while  all  of  the  parameter  files,  includ¬ 
ing  master_calib_file,dat  and  the  screening  and  parameter  files 
(SCRN_CAL  and  PARM_CAL;  §  4.2),  are  stored  in  the  di¬ 
rectory  v3.2/parmfiles. 

4.  STEP  BY  STEP 

In  this  section,  we  discuss  the  pipeline  subroutine  by  sub¬ 
routine.  Our  goal  is  to  describe  the  algorithms  employed  by 
each  subroutine  and  any  shortcomings  or  caveats  of  which  the 
user  should  be  aware. 

4.1.  OPUS 

The  Operations  Pipeline  Unified  System  (OPUS)  is  the  data- 
processing  system  used  by  the  Space  Telescope  Science  Insti¬ 
tute  to  reduce  science  data  from  the  Hubble  Space  Telescope 
( HST ).  We  use  a  /’TAL-specific  version  of  OPUS  to  manage 
our  data  processing  (Rose  et  al.  1998).  OPUS  ingests  the  data 
downlinked  by  the  spacecraft  and  produces  the  data  files  that 
serve  as  input  to  the  CalFUSE  pipeline.  OPUS  then  manages 
the  execution  of  the  pipeline  and  the  files  produced  by  CalFUSE 
and  calls  the  additional  routines  that  combine  spectra  from  each 
channel  and  exposure  into  a  set  of  observation-level  spectral 
files.  OPUS  reads  the  FUSE  Mission  Planning  Database  (which 
contains  target  information  from  the  individual  observing  pro¬ 
posals  and  instrument  configuration  and  scheduling  information 
from  the  mission  time  line)  to  populate  raw  file  header  key¬ 
words  and  to  verify  that  all  of  the  data  expected  from  an  ob¬ 
servation  were  obtained. 

OPUS  generates  six  data  files  for  each  exposure.  Four  are 
raw  data  files  (identified  by  the  suffix  "fraw.fit”),  one  for  each 
detector  segment.  One  is  a  housekeeping  file  (“hskpf.fit”)  con¬ 
taining  time-dependent  spacecraft  engineering  data.  Included 
in  this  file  are  detector  voltages,  count  rates,  and  spacecraft¬ 
pointing  information.  The  housekeeping  file  is  used  to  generate 
a  jitter  file  (“jitrf.fit”),  which  contains  information  needed  to 
correct  the  data  for  spacecraft  motion  during  an  exposure.  De¬ 
tailed  information  on  the  format  and  contents  of  each  file  is 
provided  in  the  Appendix. 

4.2.  Generate  the  Intermediate  Data  File 

The  first  task  of  the  pipeline  is  to  convert  the  raw  data  file 
into  an  intermediate  data  file  (IDF),  which  maintains  the  data 
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in  the  form  of  a  photon-event  list.  (The  format  and  contents 
of  the  IDF  are  described  in  §  A3.)  For  data  obtained  in  time- 
tag  mode,  the  module  cf_ttag_init  merely  copies  the  ar¬ 
rival  time,  X  and  Y  detector  coordinates,  and  pulse  height  of 
each  photon  event  from  the  raw  file  to  the  TIME,  XRAW, 
YRAW,  and  PHA  arrays  of  the  IDF.  A  fifth  array,  the  photon 
weight,  is  initially  set  to  unity.  Photons  whose  X  and  Y  coor¬ 
dinates  place  them  outside  of  the  active  region  of  the  detector 
are  flagged  as  described  in  §  4.3.8.  Raw  histogram  data  are 
stored  by  OPUS  as  an  image;  the  module  cf_hist_init 
converts  each  nonzero  pixel  of  that  image  into  a  single  entry 
in  the  IDF,  with  X  and  Y  equal  to  the  pixel  coordinates  (mapped 
to  their  location  on  the  unbinned  detector),  arrival  time  set  to 
the  midpoint  of  the  exposure,  and  pulse  height  set  to  20  (pos¬ 
sible  values  range  from  0  to  31).  The  arrival  time  and  pulse 
height  are  modified  later  in  the  pipeline.  The  photon  weight  is 
set  to  the  number  of  accumulated  counts  on  the  pixel,  i.e.,  the 
number  of  photons  detected  on  that  region  of  the  detector. 

The  IDF  has  two  additional  extensions.  The  first  contains 
the  good-time  intervals  (GTIs),  a  series  of  start  and  stop  times 
(in  seconds  from  the  exposure  start  time  recorded  in  the  file 
header)  computed  by  OPUS,  when  the  data  are  thought  to  be 
valid.  For  time-tag  data,  this  extension  is  copied  directly  from 
the  raw  data  file.  For  histogram  data,  a  single  GTI  is  generated 
with  START  =  0  and  STOP  =  EXPTIME  (the  exposure  time 
computed  by  OPUS).  The  final  extension  is  called  the  time¬ 
line  table  and  consists  of  16  arrays  containing  status  flags  and 
spacecraft-position,  detector  high-voltage,  and  count-rate  pa¬ 
rameters  tabulated  once  per  second  throughout  the  exposure. 
Only  the  day/night  and  OPUS  bits  of  the  time-dependent  status 
flags  are  populated  (§  A3);  the  others  are  initialized  to  zero. 
The  spacecraft-position  parameters  are  computed  as  described 
below.  The  detector  voltages  and  the  values  of  various  counters 
are  read  from  the  housekeeping  data  file. 

A  critical  step  in  the  initialization  of  the  IDF  is  populating 
the  file  header  keywords  that  describe  the  spacecraft’s  orbit 
and  control  the  subsequent  actions  of  the  pipeline.  The  names 
of  all  calibration  files  to  be  used  by  the  pipeline  are  read  from 
the  file  master_calib_file.dat  and  written  to  file  header  key¬ 
words.  (Keywords  for  each  calibration  file  are  included  in  the 
discussion  that  follows.)  Three  sets  of  calibration  files  are  time 
dependent:  the  effective  area  is  interpolated  from  the  two  files 
with  effective  dates  immediately  preceding  and  following  the 
exposure  start  time  (these  file  names  are  stored  in  the  header 
keywords  AEFF1CAL  and  AEFF2CAL);  the  scattered-light 
model  is  taken  from  the  file  with  an  effective  date  immediately 
preceding  the  exposure  start  time  (keyword  BKGD_CAL);  and 
the  orbital  elements  are  read  from  the  FUSE.TLE  file,  an  ASCII 
file  containing  NORAD  two-line  elements  for  each  day  of  the 
mission.  These  two-line  elements  are  used  to  populate  both  the 
orbital  ephemeris  keywords  in  the  IDF  header  and  the  various 
spacecraft-position  arrays  in  the  time-line  table.  Finally,  a  series 
of  data-processing  keywords  is  set  to  either  PERFORM  or 
OMIT  the  subsequent  steps  of  the  pipeline.  Once  a  step  is 


performed,  the  corresponding  keyword  is  set  to  COMPLETE. 
Some  user  control  of  the  pipeline  is  provided  by  the  screening 
and  parameter  files  (SCRN_CAL  and  PARM_CAL),  which 
allow  one,  for  example,  to  select  only  nighttime  data  or  to  turn 
off  background  subtraction.  An  annotated  list  of  file  header 
keywords,  including  the  calibration  files  used  by  the  pipeline, 
is  provided  in  the  FUSE  Instrument  and  Data  Handbook. 

Caveats:  Occasionally,  photon  arrival  times  in  raw  time-tag 
data  files  are  corrupted.  When  this  happens,  some  fraction  of 
the  photon  events  have  identical,  enormous  TIME  values,  and 
the  good-time  intervals  contain  an  entry  with  START  and 
STOP  set  to  the  same  large  value.  The  longest  valid  exposure 
spans  55  ks  (although  most  are  ~2  ks  long).  If  an  entry  in  the 
GTI  table  exceeds  this  value,  the  corresponding  entry  in  the 
time-line  table  is  flagged  as  bad  (using  the  “photon  arrival  time 
unknown”  flag;  §  A3).  Bad  TIME  values  less  than  55  ks  will 
not  be  detected  by  the  pipeline. 

Raw  histogram  files  may  also  be  corrupted.  OPUS  fills  miss¬ 
ing  pixels  in  a  histogram  image  with  the  value  21865.  The 
pipeline  sets  the  WEIGHT  of  such  pixels  to  zero  and  flags  them 
as  bad  (by  setting  the  photon’s  “fill-data  bit”;  §  A3).  Occa¬ 
sionally,  a  single  bit  in  a  histogram  image  pixel  is  flipped, 
producing  (for  high-order  bits)  a  “hot  pixel”  in  the  image.  The 
pipeline  searches  for  pixels  with  values  greater  than  8  times 
the  average  of  their  neighbors,  identifies  the  flipped  bit,  and 
resets  it. 

One  or  more  image  extensions  may  be  missing  from  a  raw 
histogram  file  (§  A2).  If  no  extensions  are  present,  the  keyword 
EXP_STAT  in  the  IDF  header  is  set  to  —1.  Exposures  with 
nonzero  values  of  EXP_STAT  are  processed  normally  by  the 
pipeline,  but  are  not  included  in  the  observation-level  spectral 
files  ultimately  delivered  to  MAST  (§  4.11).  Although  the  file 
contains  no  data,  the  header  keyword  EXPTIME  is  not  set  to 
zero. 

Early  versions  of  the  CalFUSE  pipeline  did  not  make  use 
of  the  housekeeping  files,  but  instead  employed  engineering 
information  downloaded  every  5  minutes  in  a  special  “engi¬ 
neering  snapshot”  file.  That  information  is  used  by  OPUS  to 
populate  a  variety  of  header  keywords  in  the  raw  data  file.  If 
a  housekeeping  file  is  not  available,  CalFUSE  version  3  uses 
these  keywords  to  generate  the  detector  high-voltage  and  count- 
rate  arrays  in  the  time-line  table.  Should  these  header  keywords 
be  corrupted,  the  pipeline  issues  a  warning  and  attempts  to 
estimate  the  corrupted  values.  In  such  cases,  it  is  wise  to  com¬ 
pare  the  resulting  dead-time  corrections  (§  4.3.2)  with  those  of 
other,  uncorrupted  exposures  of  the  same  target. 

4.3.  Convert  to  FARF 

The  pipeline  module  cf_convert_to_f arf  is  designed 
to  remove  detector  artifacts.  Our  goal  is  to  construct  the  data 
set  that  would  be  obtained  with  an  ideal  detector.  The  correc¬ 
tions  can  be  grouped  into  two  categories:  dead-time  effects, 
which  are  system  limitations  that  result  in  the  loss  of  photon 
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events  recorded  by  the  detector,  and  positional  inaccuracies, 
i.e.,  errors  in  the  raw  X  and  Y  pixel  coordinates  of  individual 
photon  events.  The  coordinate  system  defined  by  these  cor¬ 
rections  is  called  the  flight  alignment  reference  frame,  or  FARF. 
Corrected  coordinates  for  each  photon  event  are  written  to  the 
XFARF  and  YFARF  arrays  of  the  IDF. 

4.3.1.  Digitizer  Keywords 

The  first  subroutine  of  this  module,  cf_check_ 
digitizer,  merely  compares  a  set  of  16  IDF  file  header 
keywords  that  record  various  detector  settings  with  reference 
values  stored  in  the  calibration  file  DIGI_CAL.  Significant  dif¬ 
ferences  result  in  warning  messages  being  written  to  both  the 
file  header  and  the  exposure  trailer  file.  Such  warning  messages 
should  be  taken  seriously,  as  data  obtained  when  the  detectors 
are  not  properly  configured  are  likely  to  be  unusable.  Besides 
issuing  a  warning,  the  program  sets  the  EXP_STAT  keyword 
in  the  IDF  header  to  —2. 

4.3.2.  Detector  Dead  Time 

The  term  “dead  time”  refers  specifically  to  the  finite  time 
interval  required  by  the  detector  electronics  to  process  a  pho¬ 
ton  event.  During  this  interval,  the  detector  is  “dead”  to  in¬ 
coming  photons.  The  term  is  more  generally  applied  to  any 
loss  of  data  that  is  count-rate  dependent.  There  are  three  major 
contributions  to  the  effective  detector  dead  time  on  FUSE. 
The  first  is  due  to  limitations  in  the  detector  electronics,  which 
at  high  count  rates  may  not  be  able  to  process  photon  events 
as  fast  as  they  arrive.  The  correction  for  this  effect  is  com¬ 
puted  separately  for  each  segment  from  the  count  rate  mea¬ 
sured  at  the  detector  anode  by  the  Fast  Event  Counter  (FEC) 
and  recorded  to  the  engineering  data  stream,  typically  once 
every  16  s.  The  functional  form  of  the  correction  was  provided 
by  the  detector  development  group  at  the  University  of  Cal¬ 
ifornia,  Berkeley,  and  its  numerical  constants  were  determined 
from  in-flight  calibration  data.  It  is  applied  by  the  subroutine 
cf_electronics_dead_time. 

A  second  contribution  to  the  dead  time  comes  from  the  way 
that  the  Instrument  Data  System  processes  counts  coming  from 
the  detector.  The  IDS  can  accept  at  most  8000  counts  s  1  in 
time-tag  mode  and  32,000  counts  s  !  in  histogram  mode  from 
the  four  detector  segments  (combined).  At  higher  count  rates, 
photon  events  are  lost.  To  correct  for  such  losses,  the  subroutine 
cf_ids_dead_time  compares  the  Active  Image  Counter 
(AIC)  count  rate,  measured  at  the  back  end  of  the  detector 
electronics,  with  the  maximum  allowed  rate.  The  IDS  dead¬ 
time  correction  is  the  ratio  of  these  two  numbers  (or  unity, 
whichever  is  greater). 

A  third  contribution  occurs  when  time-tag  data  are  bundled 
into  64  kbyte  data  blocks  in  the  IDS  bulk  memory.  This 
memory  is  organized  as  a  software  FIFO  (first-in,  first-out) 
memory  buffer,  and  the  maximum  data  transfer  rate  from  it 


to  the  spacecraft  recorder  (the  FIFO  drain  rate)  is  approxi¬ 
mately  3500  events  s  \  At  higher  count  rates,  the  FIFO  will 
eventually  fill,  resulting  in  the  loss  of  one  or  more  data  blocks. 
The  effect  appears  as  a  series  of  data  dropouts,  each  a  few 
seconds  in  length,  in  the  raw  data  files.  The  correction,  com¬ 
puted  by  the  subroutine  cf_fif  o_dead_time,  is  simply  the 
ratio  of  the  AIC  count  rate  to  the  FIFO  drain  rate.  When  trig¬ 
gered,  this  correction  incorporates  (and  replaces)  the  IDS  cor¬ 
rection  discussed  above. 

The  total  dead-time  correction  (always  >1.0)  is  simply  the 
product  of  the  detector  electronics  and  IDS  corrections.  It  is 
computed  (by  the  subroutine  cf_apply_dead_time)  once 
each  second  and  applied  to  the  data  by  scaling  the  WEIGHT 
associated  with  each  photon  event.  The  mean  value  of  the 
detector  electronics,  IDS,  and  total  dead-time  corrections  are 
stored  in  the  DET_DEAD,  IDS_DEAD,  and  TOT_DEAD 
header  keywords,  respectively.  Other  possible  sources  of  dead 
time,  such  as  losses  due  to  the  finite  response  time  of  the  MCPs, 
have  a  much  smaller  effect  and  are  ignored. 

Caveats'.  Our  dead-time  correction  algorithms  are  inappro¬ 
priate  for  very  bright  targets.  If  the  header  keyword 
TOT_DEAD  >  1.5,  then  the  exposure  should  not  be  considered 
photometric.  If  DET_DEAD  >  1.5  or  if  IDS_DEAD  >1.5  and 
its  value  is  dubious  (the  housekeeping  file  for  that  exposure  is 
missing  and  the  file  header  keywords  from  which  the  count 
rates  are  estimated  appear  to  be  corrupted),  then  the  dead-time 
correction  is  assumed  to  be  meaningless  and  is  set  to  unity. 
Warning  messages  are  written  to  the  file  header  and  the  trailer 
file. 

4.3.3.  Temperature-dependent  Changes  in  Detector 
Coordinates 

The  X  and  Y  coordinates  of  a  photon  event  do  not  correspond 
to  a  physical  pixel  on  the  detector,  but  are  calculated  from 
timing  and  voltage  measurements  of  the  incoming  charge  cloud 
(Siegmund  et  al.  1997;  Sahnow  et  al.  2000a).  As  a  result,  the 
detector  coordinate  system  is  subject  to  drifts  in  the  detector 
electronics  caused  by  temperature  changes  and  other  effects. 
To  track  these  drifts,  two  signals  are  periodically  injected  into 
the  detector  electronics.  These  “stim  pulses”  appear  near  the 
upper  left  and  upper  right  corners  of  each  detector,  outside  of 
the  active  spectral  region.  The  stim  pulses  are  well  placed  for 
tracking  changes  in  the  scale  and  offset  of  the  X  coordinate, 
but  they  are  not  separated  enough  in  Y  to  track  scale  changes 
along  that  axis.  The  subroutine  cf_thermal_distort  de¬ 
termines  the  X  and  Y  centroids  of  the  stim  pulses,  computes 
the  linear  transformation  necessary  to  move  them  to  their  ref¬ 
erence  positions,  and  applies  that  transformation  to  the  X  and 
Y  coordinates  of  each  photon  event  in  the  regions  of  the  stim 
pulses  and  in  the  active  region  of  the  detector.  Events  falling 
within  64  pixels  (in  X  and  Y)  of  the  expected  stim-pulse  po¬ 
sitions  are  flagged  by  setting  the  stim-pulse  bit  in  the 
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Fig.  4. — Image  of  detector  1A  in  raw  X-  and  F-coordinates  showing  geometric  distortion.  The  image  shows  only  a  portion  of  the  detector.  It  was  constructed 
from  three  separate  exposures  with  stars  in  the  HIRS,  MDRS,  and  LWRS  apertures  of  the  LiF  1A  channel. 


LOC_FLGS  array  (§  A3).  In  raw  histogram  files,  the  stim 
pulses  are  stored  in  a  pair  of  image  extensions.  If  either  of 
these  extensions  is  missing,  the  pipeline  reads  the  expected 
positions  of  the  stim  pulses  from  the  calibration  file  STIM_CAL 
and  applies  the  corresponding  correction.  This  works  (to  first 
order)  because  the  stim  pulses  drift  slowly  with  time,  although 
short-timescale  variations  cannot  be  corrected  if  the  stim  pulses 
are  absent. 

4.3.4.  Count-Rate-dependent  Changes  in  Detector  Y  Scale 

For  reasons  not  yet  understood,  the  detector  Y  scale  varies 
with  the  count  rate,  in  the  sense  that  the  detector  image  for  a 
high-count-rate  exposure  is  expanded  in  Y.  To  measure  this 
effect,  we  tabulated  the  positions  of  individual  detector  features 
(particularly  bad-pixel  regions)  as  a  function  of  the  FEC  count 
rate  (§  4.3.2)  and  determined  the  Y  corrections  necessary  to 
shift  each  detector  feature  to  its  observed  position  in  a  low- 
count-rate  exposure.  From  this  information,  we  derived  the 
calibration  file  RATE_CAL  for  each  detector  segment.  The 
correction  is  stored  as  a  two-dimensional  image:  the  first  di¬ 
mension  represents  the  count  rate,  and  the  second  is  the  ob¬ 
served  Y  pixel  value.  The  value  of  each  image  pixel  is  the  Y 
shift  (in  pixels)  necessary  to  move  a  photon  to  its  corrected 
position.  The  subroutine  cf_count_rate_y_distort  ap¬ 
plies  this  correction  to  each  photon  event  in  the  active  region 
of  the  detector.  For  time-tag  data,  the  FEC  count  rate  is  used 
to  compute  a  time-  and  T-dependent  correction;  for  histogram 
data,  the  weighted  mean  of  the  FEC  count  rate  is  used  to  derive 
a  set  of  shifts  that  depends  only  on  Y. 

4.3.5.  Time-dependent  Changes  in  Detector  Coordinates 

As  the  detector  and  its  electronics  age,  their  properties 
change,  resulting  in  small  drifts  in  the  computed  coordinates 


of  photon  events.  These  changes  are  most  apparent  in  the  Ly/3 
airglow  features  observed  in  each  of  the  three  apertures  of  the 
LiF  and  SiC  channels  (Fig.  3),  which  drift  slowly  apart  in  Y 
as  the  mission  progresses,  indicating  a  time-dependent  stretch 
in  the  detector  Y  scale.  To  correct  for  this  stretch,  the  subroutine 
cf_time_xy_distort  applies  a  time-dependent  correction 
(stored  in  the  calibration  file  TMXY_CAL)  to  the  Y  coordinate 
of  each  photon  event  in  the  active  region  of  the  detector. 

Caveats:  Although  there  is  likely  to  be  a  similar  change  to 
the  X  coordinate,  no  measurement  of  time-dependent  drifts  in 
that  dimension  is  available,  so  no  correction  is  applied. 

4.3.6.  Geometric  Distortion 

In  an  image  of  the  detector  generated  from  raw  X  and  Y 
coordinates,  the  spectrum  is  not  straight,  but  wiggles  in  the  Y 
dimension  (Fig.  4).  To  map  these  geometric  distortions,  we 
made  use  of  two  wire  grids  (the  so-called  quantum  efficiency 
and  plasma  grids)  that  lie  in  front  of  each  detector  segment. 
Both  grids  are  regularly  spaced  and  cover  the  entire  active  area 
of  the  detectors.  Although  designed  to  be  invisible  in  the 
spectra,  they  cast  sharp  shadows  on  the  detector  when  illu¬ 
minated  directly  by  onboard  stimulation  (or  “stim”)  lamps. 
We  determined  the  shifts  necessary  to  straighten  these  shad¬ 
ows.  Their  spacing  is  approximately  1  mm,  too  great  to  mea¬ 
sure  fine-scale  structure  in  the  X  dimension,  but  sufficient  for 
the  Y  distortion.  Geometric  distortions  in  the  X  dimension 
have  the  effect  of  compressing  and  expanding  the  spectrum 
in  the  dispersion  direction,  so  the  X  distortion  correction  is 
derived  in  parallel  with  the  wavelength  calibration,  as  de¬ 
scribed  in  §  5.1.  The  geometric  distortion  corrections  are 
stored  in  a  set  of  calibration  files  (GEOM_CAL)  as  pairs  of 
16,384  x  1024  images,  one  each  for  the  X  and  Y  corrections. 
The  value  of  each  image  pixel  is  the  shift  necessary  to  move 
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a  photon  to  its  corrected  position.  This  shift  is  applied  by  the 
subroutine  cf_geometric_distort. 

Caveats:  Although  designed  to  be  invisible,  the  wire  grids 
can  cast  shadows  that  are  visible  in  the  spectra  of  astrophysical 
targets.  These  shadows  are  the  worms  discussed  in  §  6.3. 

4.3. 7.  Pulse-Height  Variations  in  Detector  X  Scale 

The  FUSE  detectors  convert  each  ultraviolet  photon  into  a 
shower  of  electrons,  for  which  the  detector  electronics  calculate 
the  X  and  Y  coordinates  and  the  total  charge,  or  pulse  height. 
Prolonged  exposure  to  photons  causes  the  detectors  to  become 
less  efficient  at  this  photon-to-electron  conversion  (a  phenom¬ 
enon  called  “gain  sag’’),  and  the  mean  pulse  height  slowly 
decreases.  Unfortunately,  the  X  coordinate  of  low-pulse-height 
photon  events  is  systematically  miscalculated  by  the  detector 
electronics.  As  the  pulse  height  decreases  with  time,  spectral 
features  appear  to  “walk”  across  the  detector.  The  strength  of 
the  effect  depends  on  the  cumulative  photon  exposure  expe¬ 
rienced  by  each  pixel  and  therefore  varies  with  location  on  the 
detector. 

To  measure  the  error  in  X  as  a  function  of  pulse  height,  we 
used  data  from  long  stim  lamp  exposures  to  construct  a  series 
of  32  detector  images,  each  containing  events  with  a  single 
pulse  height  (allowed  values  range  from  0  to  31).  We  stepped 
through  each  image  in  X,  computing  the  shift  (AX)  necessary 
to  align  the  shadow  of  each  grid  wire  with  the  corresponding 
shadow  in  a  standard  image  constructed  from  photon  events 
with  pulse  heights  between  16  and  20.  The  shifts  were 
smoothed  to  eliminate  discontinuities  and  stored  in  calibration 
files  (PHAX_CAL)  as  a  two-dimensional  image:  the  first  di¬ 
mension  represents  the  observed  X  coordinate,  and  the  second 
is  the  pulse  height.  The  value  of  each  image  pixel  is  the  walk 
correction  (AX)  to  be  added  to  the  observed  value  of  X.  This 
correction,  assumed  to  be  independent  of  detector  Y  position, 
is  applied  by  the  subroutine  cf_pha_x_distort. 

Caveats:  For  time-tag  data,  the  walk  correction  is  straight¬ 
forward  and  reasonably  accurate.  For  histogram  data,  pulse- 
height  information  is  unavailable,  so  the  subroutine 
cf_modif y_hist_pha  assigns  to  each  photon  event  the 
mean  pulse  height  for  that  aperture,  derived  from  contempo¬ 
raneous  time-tag  observations  and  stored  in  the  calibration  file 
PHAH_CAL.  While  this  trick  places  time-tag  and  histogram 
data  on  the  same  overall  wavelength  scale,  small-scale  coor¬ 
dinate  errors  due  to  localized  regions  of  gain  sag  (e.g.,  around 
bright  airglow  lines,  particularly  Ly/3 )  remain  uncorrected  in 
histogram  spectra. 

4.3.8.  Detector  Active  Region 

When  the  IDF  is  first  created,  photon  events  with  coordinates 
outside  the  active  region  of  the  detector  are  flagged  as  bad 
(§  4.2).  Once  their  coordinates  are  converted  to  the  FARE  the 
subroutine  cf_active_region  flags  as  bad  any  photons 
that  have  been  repositioned  beyond  the  active  region  of  the 


detector.  These  limits  are  read  from  the  electronics  calibration 
file  (stored  under  the  header  keyword  ELEC_CAL).  Allowed 
values  are  800  <  XFARF  <  15,583  and  0  <  YFARF  <  1023. 
The  active-area  bit  is  written  to  the  LOC_FFGS  array. 

4.3.9.  Uncorrected  Detector  Effects 

CalFUSE  does  not  perform  any  sort  of  flat-field  correction. 
Prelaunch  flat-field  observations  differ  sufficiently  from  in¬ 
flight  data  to  make  them  unsuitable  for  this  purpose,  and  in¬ 
flight  flat-field  data  are  unavailable.  (Even  if  such  data  were 
available,  any  flat-field  correction  would  be  only  approximate, 
because  MCPs  do  not  have  physical  pixels  for  which  pixel-to- 
pixel  variations  can  be  clearly  delineated;  §  4.3.3).  As  a  result, 
detector  fixed-pattern  noise  limits  the  signal-to-noise  ratio 
achievable  in  observations  of  bright  targets.  To  the  extent  that 
grating,  mirror,  and  spacecraft  motions  shift  the  spectrum  on 
the  detector  during  an  exposure,  these  fixed-pattern  features 
can  be  averaged  out.  For  some  targets,  we  deliberately  move 
the  FPAs  between  exposures  to  place  the  spectrum  on  different 
regions  of  the  detector.  Combining  the  extracted  spectra  from 
these  exposures  can  significantly  improve  the  resulting  signal- 
to-noise  ratio  (§  4.5.5).  Other  detector  effects  (including  the 
moire  pattern  discussed  in  §  6.4)  are  described  in  the  FUSE 
Instrument  and  Data  Handbook. 

4.4.  Screen  Photons 

The  module  cf_screen_photons  calls  subroutines  de¬ 
signed  to  identify  periods  of  potentially  bad  data,  such  as  Earth 
limb-angle  violations.  South  Atlantic  Anomaly  (SAA)  pas¬ 
sages,  and  detector  bursts.  A  distinct  advantage  of  CalFUSE 
version  3  over  earlier  versions  of  the  pipeline  is  that  bad  data 
are  not  discarded,  but  merely  flagged,  allowing  users  to  modify 
their  selection  criteria  without  having  to  reprocess  the  data.  To 
speed  up  processing,  the  pipeline  calculates  the  various  screen¬ 
ing  parameters  once  per  second  throughout  the  exposure,  sets 
the  corresponding  flags  in  the  STATUS_FLAGS  array  of  the 
time-line  table,  then  copies  the  flags  from  the  appropriate  entry 
of  the  time-line  table  into  the  TIMEFLGS  array  for  each  photon 
event  (§  A3).  Many  of  the  screening  parameters  applied  by  the 
pipeline  are  set  in  the  screening  parameter  file  (SCRN_CAL). 
Other  parameters  are  stored  in  various  calibration  files,  as  de¬ 
scribed  below. 

4.4.1.  Airglow  Events 

Numerous  geocoronal  emission  features  lie  within  the  FUSE 
wave  band  (Feldman  et  al.  2001).  While  the  pipeline  processes 
airglow  photons  in  the  same  manner  as  all  other  photon  events 
in  the  target  aperture,  it  is  occasionally  helpful  to  exclude  from 
consideration  regions  of  the  detector  likely  to  be  contaminated 
by  geocoronal  or  scattered  solar  emission.  These  regions  are 
listed  in  the  calibration  file  AIRG_CAL;  the  subroutine 
cf_screen_airglow  flags  as  airglow  (by  setting  the  air¬ 
glow  bit  of  the  LOC_FLGS  array  in  the  photon-event  list)  all 
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photon  events  falling  within  the  tabulated  regions — even  for 
data  obtained  during  orbital  night,  when  many  airglow  features 
are  absent. 

4.4.2.  Limb  Angle 

Spectra  obtained  when  a  target  lies  near  the  Earth’s  limb  are 
contaminated  by  scattered  light  from  strong  geocoronal  Lya 
and  O  i  emission.  To  minimize  this  effect,  the  subroutine 
cf_screen_limb_angle  reads  the  LIMB_ANGLE  array 
of  the  time-line  table,  identifies  periods  when  the  target  violates 
the  limb-angle  constraint,  and  sets  the  corresponding  flag  in 
the  STATUS_FLAGS  array  of  the  time-line  table.  Minimum 
limb  angles  for  day  and  night  observations  are  read  from  the 
BRITLIMB  and  DARKLIMB  keywords  of  the  screening  pa¬ 
rameter  file  and  copied  to  the  IDF  file  header.  The  default  limits 
are  15°  during  orbital  day  and  10°  during  orbital  night. 

4.4.3.  SAA  Passages 

The  South  Atlantic  Anomaly  marks  a  depression  in  the 
Earth’s  magnetic  field  that  allows  particles  trapped  in  the  Van 
Allen  belts  to  reach  low  altitudes.  The  high  particle  flux  in  this 
region  raises  the  background  count  rate  of  the  FUSE  detectors 
to  unacceptable  levels.  The  subroutine  cf_screen_saa 
compares  the  spacecraft’s  ground  track,  recorded  in  the  LON¬ 
GITUDE  and  LATITUDE  arrays  of  the  time-line  table,  with 
the  limits  of  the  SAA  (stored  in  the  calibration  file  SAAC_CAL 
as  a  binary  table  of  latitude-longitude  pairs)  and  flags  as  bad 
periods  when  data  were  taken  within  the  SAA.  Our  SAA  model 
was  derived  from  orbital  information  and  onboard  counter  data 
from  the  first  3  years  of  the  FUSE  mission. 

Caveats:  Because  the  SAA  particle  flux  is  great  enough  to 
damage  the  FUSE  detectors,  we  end  most  exposures  before 
entering  the  SAA  and  lower  the  detector  high  voltage  during 
each  SAA  pass.  As  a  result,  very  little  data  is  actually  flagged 
by  this  screening  step. 

4.4.4.  Detector  High  Voltage 

The  detector  high  voltage  is  set  independently  for  each 
detector  segment  (1A,  IB,  2A,  2B).  During  normal  operations, 
the  voltage  on  each  segment  alternates  between  its  nominal 
full  voltage  and  a  reduced  SAA  level.  The  SAA  level  is  low 
enough  that  the  detectors  are  not  damaged  by  the  high  count 
rates  that  result  from  SAA  passes,  and  it  is  often  used  between 
science  exposures  to  minimize  detector  exposure  to  bright 
airglow  emission.  The  full-voltage  level  is  the  normal  oper¬ 
ating  voltage  used  during  science  exposures.  It  is  raised  reg¬ 
ularly  to  compensate  for  the  effects  of  detector  gain  sag. 
Without  this  compensation,  the  mean  pulse  height  of  real 
photon  events  would  gradually  fall  below  our  detection 
threshold.  Unfortunately,  there  is  a  limit  above  which  the  full- 
voltage  level  cannot  be  raised.  Detector  segment  2A  reached 
this  limit  in  2003  and  has  not  been  raised  since;  it  will  grad¬ 
ually  become  less  sensitive  as  the  fraction  of  low-pulse-height 


events  increases.  The  subroutine  cf_screen_high_ 
voltage  reads  the  instantaneous  value  of  the  detector  high 
voltage  from  the  HIGH_VOLTAGE  array  of  the  time-line  table, 
compares  it  with  the  nominal  full-voltage  level  (stored  as  a 
function  of  time  in  the  calibration  file  VOLT_CAL),  and  flags 
periods  of  low  voltage  as  bad. 

For  any  number  of  reasons,  an  exposure  can  be  obtained 
with  the  detector  high  voltage  at  less  than  the  full-voltage  level. 
To  preserve  as  much  of  this  data  as  possible,  we  examined  all 
of  the  low-voltage  exposures  taken  during  the  first  4  years  of 
the  mission  and  found  that  for  detector  segments  1A,  IB,  and 
2B,  the  data  quality  is  good  whenever  the  detector  high  voltage 
is  greater  than  85%  of  the  nominal  (time  dependent)  full-volt- 
age  level.  For  segment  2A,  data  obtained  with  the  high  voltage 
greater  than  90%  of  full  are  good,  lower  than  80%  are  bad, 
and  between  80%  and  90%  are  of  variable  quality.  In  this 
regime,  the  pipeline  flags  the  affected  data  as  good  but  writes 
warning  messages  to  both  the  IDF  header  and  the  trailer  file. 
When  this  warning  is  present  in  time-tag  data,  the  user  should 
examine  the  distribution  of  pulse  heights  in  the  target  aperture 
to  ensure  that  the  photon  events  are  well  separated  from  the 
background  (§  4.4.12).  For  histogram  data,  the  spectral  pho¬ 
tometry  and  wavelength  scale  are  most  likely  to  be  affected. 

Caveats:  If  the  header  keywords  indicate  that  the  detector 
voltage  was  high,  low,  or  changed  during  an  exposure,  the  IDF 
initialization  routines  (§  4.2)  write  a  warning  message  to  the 
trailer  file.  If  a  valid  housekeeping  file  is  available  for  the 
exposure,  this  warning  may  be  safely  ignored,  because  the 
pipeline  uses  housekeeping  information  to  populate  the 
HIGH_VOLTAGE  array  in  the  time-line  table  and  properly 
excludes  time  intervals  when  the  voltage  was  low.  If  the  house¬ 
keeping  file  is  not  present,  each  entry  of  the  HIGH_VOLTAGE 
array  is  set  to  the  “HV  bias  maximum  setting”  reported  in  the 
IDF  header.  In  this  case,  the  pipeline  has  no  information  about 
time-dependent  changes  in  the  detector  high  voltage,  and  warn¬ 
ings  about  voltage  level  changes  should  be  investigated  by  the 
user. 

4.4.5.  Event  Bursts 

Occasionally,  the  FUSE  detectors  register  large  count  rates 
for  short  periods  of  time.  These  event  bursts  can  occur  on  one 
or  more  detectors  and  often  have  a  complex  distribution  across 
the  detector,  including  scalloping  and  sharp  edges  (Fig.  5). 
CalFUSE  includes  a  module  that  screens  the  data  to  identify 
and  exclude  bursts.  The  subroutine  cf_screen_burst  com¬ 
putes  the  time-dependent  count  rate  using  data  from  back¬ 
ground  regions  of  the  detector  (excluding  airglow  features)  and 
applies  a  median  filter  to  reject  time  intervals  whose  count  rates 
differ  by  more  than  5  standard  deviations  (the  value  can  be  set 
by  the  user)  from  the  mean.  The  algorithm  rejects  any  time 
interval  in  which  the  background  rate  rises  rapidly,  as  when 
an  exposure  extends  into  an  SAA  or  the  target  nears  the  Earth 
limb.  The  background  rate  computed  by  the  burst-rejection 


2007  PASP,  119:527-555 


536  DIXON  ET  AL. 


Fig.  5. — Segments  of  detector  1A  images  showing  filamentary  {top)  and  checkerboard  {bottom)  bursts.  Checkerboard  bursts  typically  fill  the  detector,  save  for 
the  region  around  the  LiF  Ly/3  lines  on  detector  1A. 


algorithm  is  stored  in  the  BKGD_CNT_RATE  array  of  the 
time-line  table  and  included  on  the  count-rate  plots  generated 
for  each  exposure  (§  4.10).  Burst  rejection  is  possible  only  for 
data  obtained  in  time-tag  mode. 

Caveats:  Careful  examination  of  long  background  obser¬ 
vations  reveals  that  many  are  contaminated  by  emission  from 
bursts  too  faint  to  trigger  the  burst-detection  algorithm.  Ob¬ 
servers  studying,  for  example,  diffuse  emission  from  the  in¬ 
terstellar  medium  should  be  alert  to  the  possibility  of  such 
contamination. 

4.4.6.  Spacecraft  Drift 

Pointing  of  the  FUSE  spacecraft  was  originally  controlled 
with  four  reaction  wheels,  which  typically  maintained  a  point¬ 
ing  accuracy  of  0.2"-0.3".  In  late  2001,  two  of  the  reaction 
wheels  failed,  and  it  became  necessary  to  control  the  spacecraft 
orientation  along  one  axis  with  magnetic  torquer  bars.  The 
torquer  bars  can  exert  only  about  10%  of  the  force  produced 
by  the  reaction  wheels,  and  the  available  force  depends  on  the 
strength  and  relative  orientation  of  the  Earth’s  magnetic  field. 
Thus,  spacecraft  drift  increased  dramatically  along  this  axis, 
termed  the  antisymmetric  or  A  axis.  Drifts  about  the  A  axis 
shift  the  spectra  of  point-source  targets  in  a  direction  45°  from 
the  dispersion  direction  (i.e.,  AX  =  AF).  These  motions  can 
substantially  degrade  the  resolution  of  the  spectra,  so  proce¬ 


dures  have  been  implemented  to  correct  the  data  for  spacecraft 
motion  during  an  exposure.  For  time-tag  observations  of  point 
sources,  we  reposition  individual  photon  events.  For  histogram 
observations,  we  correct  only  for  the  exposure  time  lost  to  large 
excursions  of  the  spacecraft.  The  ability  to  correct  for  spacecraft 
drift  became  even  more  important  when  a  third  reaction  wheel 
failed  in  2004  December. 

The  correction  of  photon-event  coordinates  for  spacecraft 
motion  is  discussed  in  §  4.5.7.  During  screening,  the  subroutine 
cf_screen_j  itter  merely  flags  times  when  the  target  is 
out  of  the  aperture,  defined  as  those  for  which  either  AX  or 
AF,  the  pointing  error  in  the  dispersion  or  cross-dispersion 
direction,  respectively,  is  greater  than  30",  the  width  of  the 
LWRS  aperture.  These  limits,  set  by  the  keywords  DX_MAX 
and  DY_MAX  in  the  CalFUSE  parameter  file  (PARM_CAL), 
underestimate  the  time  lost  to  pointing  excursions,  but  smaller 
limits  can  lead  to  the  rejection  of  good  data  for  some  channels. 
Also  flagged  as  bad  are  times  when  the  jitter-tracking  flag 
TRKFLG  =  —1,  indicating  that  the  spacecraft  is  not  tracking 
properly.  If  TRKFLG  =  0,  no  tracking  information  is  available 
and  no  times  are  flagged  as  bad.  Pointing  information  is  read 
from  the  jitter  file  (JITR_CAL;  §  A2).  If  the  jitter  file  is  not 
present  or  the  header  keyword  JIT_STAT  =  1  (indicating  that 
the  jitter  file  is  corrupted),  cf_screen_j  itter  issues  a 
warning  and  exits;  again,  no  times  are  flagged  as  bad. 
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4.4.7.  User-defined  Good-Time  Intervals 

One  bit  of  the  status  array  is  reserved  for  user-defined  GTIs. 
For  example,  to  extract  data  corresponding  to  a  particular  phase 
of  a  binary  star  orbit,  one  would  flag  data  from  all  other 
phases  as  bad.  A  number  of  tools  exist  to  set  this  flag,  in¬ 
cluding  cf_ed.it  (available  from  MAST).  CalFUSE  users 
can  specify  good-time  intervals  by  setting  the  appropriate 
keywords  (NUSERGTI,  GTIBEG01,  GTIEND01,  etc.)  in  the 
screening  parameter  file.  (Times  are  in  seconds  from  the  ex¬ 
posure  start  time.)  If  these  keywords  are  set,  the  subroutine 
cf_set_user_gtis  flags  times  outside  of  these  good-time 
intervals  as  bad. 

4.4.8.  Time-dependent  Status  Flags 

Once  the  status  flags  in  the  time-line  table  are  populated, 
the  subroutine  cf_set_photon_flags  copies  them  to  the 
corresponding  entries  in  the  photon-event  list.  For  time-tag 
data,  this  process  is  straightforward:  match  the  times  and  copy 
the  flags.  Fleader  keywords  in  the  IDF  record  the  number  of 
photon  events  falling  in  bad-time  intervals  or  outside  of  the 
detector  active  area;  the  number  of  seconds  lost  to  bursts,  SAAs, 
etc.;  and  the  remaining  night  exposure  time.  If  more  than  90% 
of  the  exposure  is  lost  to  a  single  cause,  an  explanatory  note 
is  written  to  the  trailer  file. 

The  task  is  more  difficult  for  histogram  data,  for  which  pho¬ 
ton-arrival  information  is  unavailable.  We  distinguish  between 
time  flags  that  represent  periods  of  lost  exposure  time  (low 
detector  voltage  or  target  out  of  aperture)  and  those  that  rep¬ 
resent  periods  of  data  contamination  (limb-angle  violations  or 
SAAs).  For  the  former,  we  need  only  modify  the  exposure  time; 
for  the  latter,  we  must  flag  the  exposure  as  being  contaminated. 
Our  goal  is  to  set  the  individual  photon  flags  and  header  key¬ 
words  so  that  the  pipeline  behaves  in  the  following  way:  when 
processing  a  single  exposure,  it  treats  all  photon  events  as  good. 
When  combining  data  from  multiple  exposures,  it  excludes 
contaminated  exposures  (defined  below).  To  this  end,  we  gen¬ 
erate  an  8  bit  status  word  containing  only  day/night  informa¬ 
tion:  if  the  exposure  is  more  than  10%  day,  the  day  bit  is  set. 
This  status  word  is  copied  onto  the  time-dependent  status  flag 
of  each  photon  event.  We  generate  a  second  8  bit  status  word 
containing  information  about  limb-angle  violations  and  SAAs: 
if  a  single  second  is  lost  to  one  of  these  events,  the  corre¬ 
sponding  flag  is  set  and  a  message  is  written  to  the  trailer  file. 
(To  avoid  rejecting  an  exposure  that,  for  example,  abuts  an 
SAA,  we  ignore  its  initial  and  final  20  s  in  this  analysis.)  The 
status  word  is  stored  in  the  file  header  keyword  EXP_STAT 
(unless  that  keyword  has  already  been  set;  see  §§  4.2  and  4.3.1 ). 
The  routines  used  by  the  pipeline  to  combine  data  from  multiple 
exposures  into  a  single  spectrum  (§  4.11)  reject  data  files  in 
which  this  keyword  is  nonzero.  The  number  of  bad  events,  the 
exposure  time  lost  to  periods  of  low  detector  voltage  or  space¬ 
craft  jitter,  and  the  exposure  time  during  orbital  night  are  written 
to  the  file  header,  just  as  for  time-tag  data. 


Only  in  this  subroutine  is  the  DAYNIGHT  keyword  read 
from  the  screening  parameter  file  and  written  to  the  IDF  file 
header.  Allowed  values  are  DAY,  NIGHT,  and  BOTH.  The 
default  is  BOTH.  For  most  flags,  if  the  bit  is  set  to  1,  the  photon 
event  is  rejected.  The  day/night  flag  is  different:  it  is  always 
1  for  day  and  0  for  night.  The  pipeline  must  read  and  interpret 
the  DAYNIGHT  keyword  before  accepting  or  rejecting  an 
event  based  on  the  value  of  its  day/night  flag. 

4.4.9.  Good-Time  Intervals 

Once  the  time-dependent  screening  is  complete,  the  subrou¬ 
tine  cf_set_good_time_intervals  calculates  a  new  set 
of  good-time  intervals  from  information  in  the  time-line  table 
and  writes  them  to  the  second  extension  of  the  IDF  (§  4.2). 
For  time-tag  data,  all  of  the  TIMEFLGS  bits  are  used  and  the 
DAYNIGHT  filter  is  applied.  For  histogram  data,  the  bits  cor¬ 
responding  to  limb-angle  violations  and  SAAs  are  ignored, 
since  data  arriving  during  these  events  cannot  be  excluded.  The 
DAYNIGHT  filter  is  applied  (assuming  that  all  are  day  photons 
if  the  exposure  is  more  than  10%  day).  The  exposure  time, 
EXPTIME  =  E  (STOP  —  START),  summed  over  all  entries 
in  the  GTI  table,  is  then  written  to  the  IDF  file  header. 

4.4.10.  Histogram  Arrival  Times 

For  histogram  data,  all  of  the  photon  events  in  an  IDF  are 
initially  assigned  an  arrival  time  equal  to  the  midpoint  of  the 
exposure.  Should  this  instant  fall  in  a  bad-time  interval,  the 
data  may  be  rejected  by  a  subsequent  step  of  the  pipeline  or 
one  of  our  postprocessing  tools.  To  avoid  this  possibility,  the 
subroutine  cf_modif y_hist_times  resets  all  photon-ar¬ 
rival  times  to  the  midpoint  of  the  exposure’s  longest  good-time 
interval.  This  subroutine  is  not  called  for  time-tag  data. 

4.4.11.  Bad-Pixel  Regions 

Images  of  the  FUSE  detectors  reveal  a  number  of  dead  spots 
that  may  be  surrounded  by  a  bright  ring  (see  the  FUSE  In¬ 
strument  and  Data  Handbook  for  examples).  The  subroutine 
cf_screen_bad_pixels  reads  a  list  of  bad-pixel  regions 
from  a  calibration  file  (QUAL_CAL)  and  flags  as  bad  all  photon 
events  whose  XFARF  and  YFARF  coordinates  fall  within  the 
tabulated  limits.  A  bad-pixel  map,  constructed  later  in  the  pipe¬ 
line  (§  4.8),  is  used  by  the  optimal-extraction  algorithm  to 
correct  for  flux  lost  to  dead  spots. 

4.4.12.  Pulse-Height  Limits 

For  time-tag  data,  the  pulse  height  of  each  photon  event  is 
recorded  in  the  IDF.  Values  range  from  0  to  31  in  arbitrary 
units.  A  typical  pulse-height  distribution  has  a  peak  at  low 
values  due  to  the  intrinsic  detector  background,  a  Gaussian- 
like  peak  near  the  middle  of  the  range  due  to  “real”  photons, 
and  a  tail  of  high-pulse-height  events,  which  likely  represent 
the  superposition  of  two  photons  and  therefore  are  not  reliable. 
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In  addition,  the  detector  electronics  selectively  discard  high 
pulse-height  events  near  the  top  and  bottom  of  the  detectors 
(i.e.,  with  large  or  small  values  of  Y).  We  can  thus  improve 
the  signal-to-noise  ratio  of  faint  targets  by  rejecting  photon 
events  with  extreme  pulse-height  values.  Pulse-height  limits 
(roughly  2-24)  are  defined  for  each  detector  segment  in  the 
screening  parameter  file  (SCRN_CAL).  The  subroutine 
cf_screen_pulse_height  flags  photon  events  with 
pulse  heights  outside  of  this  range  (by  setting  the  appropriate 
bit  in  the  LOC_FLGS  array;  §  A3)  and  writes  the  pulse-height 
limits  used  and  the  number  of  photon  events  rejected  to  the 
IDF  file  header.  This  procedure  is  not  performed  on  histogram 
data. 

Caveats:  We  do  not  recommend  the  use  of  narrow  pulse- 
height  ranges  to  reduce  the  detector  background  in  FUSE  data. 
Careful  analysis  has  shown  that  limits  that  are  more  stringent 
than  the  default  values  can  result  in  significant  flux  losses  across 
small  regions  of  the  detector,  particularly  in  the  LiF  IB  channel, 
resulting  in  apparent  absorption  features  that  are  not  real. 

4.5.  Remove  Motions 

Having  corrected  the  data  for  various  detector  effects  and 
identified  periods  of  bad  data,  we  continue  to  work  backward 
through  the  instrument,  correcting  for  spectral  motions  on  the 
detector  due  to  the  movements  of  various  optical  components — 
and  even  of  the  spacecraft  itself.  This  task  is  performed  by  the 
module  cf_remove_motions.  It  begins  by  reading  the 
XFARF  and  YFARF  coordinates  of  each  photon  event  from 
the  IDF.  It  concludes  by  writing  the  motion-corrected  coor¬ 
dinates  to  the  X  and  Y  arrays  of  the  same  file, 

4.5.1.  Locate  Spectra  on  the  Detector 

The  LiF  and  SiC  channels  each  produce  three  spectra,  one 
from  each  aperture,  for  a  total  of  six  spectra  per  detector  seg¬ 
ment  (Fig.  3).  Because  motions  of  the  optical  components  can 
shift  these  spectra  on  the  detector,  the  first  step  is  to  determine 
the  Y  centroid  of  each.  To  do  this,  we  use  the  following  al¬ 
gorithm:  First,  we  project  the  airglow  photons  onto  the  Y  axis 
(summing  all  values  of  X  for  each  value  of  F)  and  search  the 
resulting  histogram  for  peaks  within  70  pixels  of  the  expected 
Y  position  of  the  LWRS  spectrum.  If  the  airglow  feature  is 
sufficiently  bright  (33  counts  in  141  Y  pixels),  we  adopt  its 
centroid  as  the  airglow  centroid  for  the  LWRS  aperture  and 
compute  its  offset  from  the  expected  value  stored  in  the 
CHID_CAL  calibration  file.  If  the  airglow  feature  is  too  faint, 
we  adopt  the  expected  centroid  and  assume  an  offset  of  zero. 
We  apply  the  offset  to  the  expected  centroids  of  the  MDRS 
and  HIRS  apertures  to  obtain  their  airglow  centroids.  Second, 
we  project  the  nonairglow  photons  onto  the  Y  axis  and  subtract 
a  uniform  background.  Airglow  features  fill  the  aperture,  but 
point-source  spectra  are  considerably  narrower  in  Y  and  need 
not  be  centered  in  the  aperture.  For  each  aperture,  we  search 
for  a  5  a  peak  within  40  pixels  of  the  airglow  centroid.  If  we 


find  it,  we  use  its  centroid;  otherwise,  we  use  the  airglow 
centroid.  This  scheme,  implemented  in  the  subroutine 
cf_find_spectra,  allows  for  the  possibility  that  an  astro- 
physical  spectrum  may  appear  in  a  nontarget  aperture. 

For  each  of  the  six  spectra,  two  keywords  are  written  to  the 
IDF  file  header:  YCENT  contains  the  computed  Y  centroid, 
and  YQUAL  contains  a  quality  flag.  The  flag  is  HIGH  if  the 
centroid  was  computed  from  a  point-source  spectrum,  ME¬ 
DIUM  if  computed  from  an  airglow  spectrum,  and  LOW  if  the 
tabulated  centroid  was  used.  It  is  possible  for  the  user  to  specify 
the  target  centroid  by  setting  the  SPEX_SIC  and  SPEC_LIF 
keywords  in  the  CalFUSE  parameter  file  (PARM_CAL).  Two 
other  keywords,  EMAX_SIC  and  EMAX_LIF,  limit  the  offset 
between  the  expected  and  calculated  centroids:  if  the  calculated 
centroid  differs  from  the  predicted  value  by  more  than  this 
limit,  the  pipeline  uses  the  default  centroid. 

Caveats:  On  detector  1,  the  SiC  LWRS  spectrum  falls  near 
the  bottom  edge  of  the  detector  (Fig.  3).  Because  the  back¬ 
ground  level  rises  steeply  near  this  edge,  the  calculated  centroid 
can  be  pulled  (incorrectly)  to  lower  values  of  Y.  especially  for 
faint  targets. 

4.5.2.  Assign  Photons  to  Channels 

The  subroutine  cf_identif y_channel  assigns  each 
photon  to  a  channel,  where  “channel”  now  refers  to  one  of  the 
six  spectra  on  each  detector  (Fig.  3).  For  each  channel,  ex¬ 
traction  windows  for  both  point-source  and  extended  targets 
are  tabulated  in  the  calibration  file  CHID_CAL,  along  with 
their  corresponding  spectral  Y  centroids.  These  extraction  limits 
encompass  at  least  99.5%  of  the  target  flux.  For  the  target 
channels,  identified  in  the  APERTURE  header  keyword,  we 
use  either  the  point-source  or  extended  extraction  windows,  as 
indicated  by  the  SRC_TYPE  keyword;  for  the  other  (presum¬ 
ably  airglow)  channels,  we  use  the  extended  extraction  win¬ 
dows.  The  offset  between  the  calculated  and  tabulated  spectral 
Y  centroids  (§  4.5.1)  is  used  to  shift  each  extraction  window 
to  match  the  data. 

To  ensure  that,  in  the  event  that  two  extraction  windows 
overlap,  photon  events  falling  in  the  overlap  region  are  assigned 
to  the  more  likely  channel,  photon  coordinates  (XFARF  and 
YFARF)  are  compared  with  the  extraction  limits  of  the  six 
spectral  channels  in  the  following  order:  first  the  target  channels 
(LiF  and  SiC),  then  the  airglow  channels  (LiF  and  SiC)  cor¬ 
responding  to  the  larger  nontarget  aperture,  and  finally  the  air¬ 
glow  channels  (LiF  and  SiC)  corresponding  to  the  smaller  non¬ 
target  aperture.  For  example,  if  the  target  were  in  the  MDRS 
aperture,  the  search  order  would  be  MDRS  LiF,  MDRS  SiC, 
LWRS  LiF,  LWRS  SiC,  HIRS  LiF,  and  HIRS  SiC.  The  process 
stops  when  a  match  is  made.  The  channel  assignment  of  each 
photon  event  is  stored  in  the  CHANNEL  array  (§  A3);  photon 
events  that  do  not  fall  in  an  extraction  window  are  assigned  a 
CHANNEL  value  of  0. 

Channel  assignment  is  performed  twice,  once  before  the 
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motion  corrections  and  once  after.  The  first  time,  all  extraction 
windows  are  padded  by  ±  1 0  Y  pixels  to  accommodate  errors 
in  the  channel  centroids.  The  second  time,  no  padding  is 
applied  to  time-tag  data.  Histogram  data,  which  are  generally 
binned  by  8  pixels  in  Y,  present  a  special  challenge:  the  geo¬ 
metric  correction  (§  4.3.6)  can  move  a  row  of  data  out  of  the 
extraction  window,  producing  a  significant  loss  of  flux.  To 
prevent  this,  histogram  extraction  windows  are  padded  by  an 
additional  ±8  f  pixels  (or  an  amount  equal  to  the  binning 
factor  in  Y,  if  other  than  8). 

4.5.3.  Track  Y  Centroids  with  Time 

For  the  LiF  and  SiC  target  apertures,  the  subroutine 
cf_calculate_ycent_motion  computes  the  spectral  Y 
centroid  as  a  function  of  time  throughout  the  exposure.  The 
algorithm  requires  500  photon  events  to  compute  an  average, 
so  the  centroid  is  updated  more  often  for  bright  targets  than 
for  faint  ones.  Photon  events  flagged  as  airglow  are  ignored. 
The  results  are  stored  in  the  YCENT_LIF  and  YCENT_SIC 
arrays  of  the  time-line  table,  but  are  not  currently  used  by  the 
pipeline.  This  calculation  is  not  performed  for  data  obtained 
in  histogram  mode. 

4.5.4.  Correct  for  Grating  Motion 

The  FUSE  spectrograph  gratings  are  subject  to  small,  ther¬ 
mally  induced  motions  on  orbital,  diurnal,  and  precessional 
(60  day)  timescales;  an  additional  long-term,  nonperiodic  drift 
is  also  apparent.  These  motions  can  shift  the  target  and  air- 
glow  spectra  by  as  much  as  15  pixels  (peak  to  peak)  in  both 
the  X  and  Y  dimensions.  Measurements  of  the  Ly/3  airglow 
line  in  thousands  of  exposures  obtained  throughout  the  mis¬ 
sion  reveal  that  the  gratings’  orbital  motion  depends  on  three 
parameters:  beta  angle  (the  angle  between  the  target  and  the 
anti-Sun  vector),  pole  angle  (the  angle  between  the  target  and 
the  orbit  pole),  and  spacecraft  roll  angle  (east  of  north,  stored 
in  the  file  header  keyword  APER_PA).  The  subroutine 
cf_grating_motion  compares  the  beta,  pole,  and  roll  an¬ 
gles  of  the  spacecraft  with  a  grid  of  values  in  the  calibration 
file  GRAT_CAL,  reads  the  appropriate  correction,  and  com¬ 
putes  the  X  and  Y  photon  shifts.  The  grating-motion  correction 
is  applied  to  all  photon  events  with  CHANNEL  >  0;  photon 
events  not  assigned  to  a  channel  are  not  moved. 

Caveats:  Some  combinations  of  beta  and  pole  angle  are  too 
poorly  sampled  for  us  to  derive  a  grating-motion  correction; 
for  these  regions,  no  correction  is  applied.  At  present,  only 
corrections  for  the  orbital  and  long-term  grating  motions  are 
available.  Because  all  photon  events  in  histogram  data  are  as¬ 
signed  the  same  arrival  time  (the  midpoint  of  the  longest  good¬ 
time  interval),  they  receive  the  same  grating-motion  correction. 

4.5.5.  Correct  for  FPA  Motion 

The  four  focal  plane  assemblies  (shown  in  Fig.  1)  can  be 
moved  independently  in  either  the  X  or  Z  direction.  FPA 


motions  in  the  X  direction  are  used  to  correct  for  mirror 
misalignments  and  to  perform  FP  splits  (described  below). 
FPA  motions  in  the  Z  direction  are  used  to  place  the  ap¬ 
ertures  in  the  focal  plane  of  the  spectrograph.  (Strictly  speak¬ 
ing,  an  FPA  moves  along  the  tangent  to  or  the  radius  of  the 
spectrograph  Rowland  circle,  not  the  X  and  Z  axes  shown 
in  Fig.  1.)  Both  motions  change  the  spectrograph  entrance 
angle,  shifting  the  target  spectrum  on  the  detector.  The  FUSE 
wavelength  calibration  is  derived  from  a  single  stellar  ob¬ 
servation  obtained  at  a  particular  FPA  position.  The  sub¬ 
routine  cf_fpa_position  computes  the  shift  in  pixels 
(AY)  necessary  to  move  each  spectrum  from  its  observed  X 
position  on  the  detector  to  that  of  the  wavelength-calibration 
target.  The  X  and  Z  positions  of  the  LiF  and  SiC  FPAs  are 
stored  in  file  header  keywords,  various  spectrograph  parameters 
are  stored  in  a  calibration  file  (SPEC_CAL),  and  the  wave¬ 
length  calibration  and  the  FPA  position  of  the  wavelength- 
calibration  target  are  stored  in  the  WAVE_CAL  file.  Shifts  are 
computed  for  both  the  LiF  and  SiC  channels.  The  appropriate 
shift  is  applied  to  all  photon  events  with  CHANNEL  >  0; 
photon  events  not  assigned  to  a  channel  are  ignored. 

The  FUSE  detectors  suffer  from  fixed-pattern  noise.  Astig¬ 
matism  in  the  instrument  spreads  a  typical  resolution  element 
over  several  hundred  detector  pixels  (predominantly  in  the 
cross-dispersion  dimension),  mitigating  this  effect,  but  to 
achieve  a  signal-to-noise  ratio  greater  than  ~30,  one  must  re¬ 
move  the  remaining  fixed-pattern  noise.  A  useful  technique  is 
the  focal-plane  split.  FP  splits  are  performed  by  obtaining  a 
series  of  MDRS  or  HIRS  exposures  at  several  FPA  X  positions. 
Moving  the  FPA  in  the  X  dimension  (and  moving  the  satellite 
to  center  the  target  in  the  aperture)  between  exposures  shifts 
both  target  and  airglow  spectra  in  the  dispersion  direction  on 
the  detector.  CalFUSE  shifts  each  spectrum  to  the  standard  X 
position  expected  by  our  wavelength  calibration  routines.  If  the 
signal-to-noise  ratio  in  the  spectra  obtained  at  each  FPA  po¬ 
sition  is  high  enough,  it  is  possible  for  the  user  to  disentangle 
the  source  spectrum  from  the  detector  fixed-pattern  noise;  how¬ 
ever,  simply  combining  extracted  spectra  obtained  at  different 
FPA  positions  will  average  out  most  of  the  small-scale  detector 
features. 

4.5.6.  Correct  for  Mirror  Motion 

The  spectrograph  mirrors  are  subject  to  thermal  motions  that 
shift  the  target’s  image  within  the  spectrograph  aperture  and 
thus  its  spectrum  in  both  X  and  Y  on  the  detector.  A  source  in 
either  of  the  SiC  channels  may  move  as  much  as  6"  in  a  period 
of  2  ks.  This  motion  has  two  effects  on  the  data:  first,  flux  will 
be  lost  if  the  source  drifts  (partially  or  completely)  out  of  the 
aperture;  second,  spectral  resolution  will  be  degraded  (for 
LWRS  observations)  as  the  spectrum  shifts  on  the  detector. 
Diffuse  sources,  such  as  airglow  emission,  fill  the  aperture,  so 
their  spectra  are  unaffected  by  mirror  motion. 

When  the  LiF  1  channel  is  used  for  guiding,  motions  of  the 
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LiF  1  mirror  are  corrected  by  the  spacecraft  itself.  Only  the 
LiF  2  and  SiC  spectra  must  be  corrected  by  CalFUSE.  In  theory, 
the  switch  from  LiF  1  to  LiF  2  as  the  primary  channel  for 
guiding  the  spacecraft  (§  2)  should  require  another  set  of  cal¬ 
ibration  files.  In  practice,  the  LiF  2  mirror  motion  in  the  dis¬ 
persion  direction  tracks  that  of  the  LiF  1  mirror.  The  mirror- 
motion  correction  is  stored  as  a  function  of  time  since  orbital 
sunset  (via  the  TIME_SUNSET  array  in  the  time-line  table)  in 
the  calibration  file  MIRR_CAL.  The  correction  (AX)  is  applied 
by  the  subroutine  cf_mirror_motion  to  all  photon  events 
within  the  target  aperture;  photon  events  in  other  apertures  and 
those  not  assigned  to  a  channel  are  ignored.  This  correction  is 
not  applied  to  extended  sources.  Because  all  photon  events  in 
histogram  data  are  assigned  the  same  arrival  time  (generally 
the  midpoint  of  the  exposure),  they  receive  the  same  mirror- 
motion  correction. 

Caveats :  We  correct  only  the  relative  mirror  motion  during 
an  orbit,  not  the  absolute  mirror  offset  based  on  longer  term 
trends.  We  do  not  correct  for  mirror  motions  in  the  Y  dimension. 
Finally,  because  the  shifts  for  the  SiC  1  and  SiC  2  mirrors  are 
similar,  we  adopt  a  single  correction  for  both  channels. 

4.5.7.  Correct  for  Spacecraft  Motion 

Spacecraft  motions  during  an  exposure  shift  the  target  spec¬ 
trum  on  the  detector  and  thus  degrade  spectral  resolution.  The 
subroutine  cf_satellite_j  itter  uses  pointing  infor¬ 
mation  stored  in  the  jitter  file  (JITR_CAL;  §  A2)  to  correct 
the  observed  coordinates  of  photon  events  for  these  motions. 
Pointing  errors  in  arcseconds  are  converted  to  X  and  Y  pixel 
shifts  and  applied  to  all  photon  events  within  the  target  ap¬ 
erture;  events  in  other  apertures  and  those  not  assigned  to  a 
channel  are  ignored.  The  correction  is  applied  only  if  the  jitter¬ 
tracking  flag  TRKFLG  >  0,  indicating  that  valid  tracking 
information  is  available.  TRKFLG  values  rise  from  1  to  5  as 
the  reliability  of  the  pointing  information  increases.  The  min¬ 
imum  acceptable  value  of  the  TRKFLG  can  be  adjusted  by 
modifying  the  TRKFLG  keyword  in  the  CalFUSE  parameter 
file  (PARM_CAL). 

4.5.8.  Recompute  Spectral  Centroids 

Once  all  spectral  motions  are  removed  from  the  data,  the 
subroutine  cf_calculate_y_centroid  recomputes  the 
spectral  Y  centroids.  Separate  source  and  airglow  centroids  are 
determined  for  each  aperture  in  turn,  from  largest  to  smallest. 
(The  former  is  meaningless  if  the  aperture  does  not  contain  a 
source.)  The  offset  between  the  measured  airglow  centroid  in 
the  LWRS  aperture  and  the  tabulated  centroid  (from  the  cali¬ 
bration  file  CHID_CAL)  is  used  to  compute  the  airglow  cen¬ 
troids  for  the  MDRS  and  HIRS  apertures;  the  computed  MDRS 
and  HIRS  airglow  centroids  are  ignored.  The  YCENT  value 
written  to  the  IDF  file  header  is  determined  by  the  quality  flag 
previously  set  by  cf_find_spectra  (§  4.5.1):  if  YQUAL  = 
HIGH,  the  source  centroid  is  used;  if  YQUAL  =  MEDIUM, 


the  airglow  centroid  is  used;  and  if  YQUAL  =  LOW,  the 
tabulated  centroid  is  used.  The  SPEX_SIC,  SPEC_LIF, 
EMAX_SIC,  and  EMAX_LIF  keywords  in  the  CalFUSE  pa¬ 
rameter  file  (PARM_CAL)  have  the  effects  discussed  in  §  4.5.1. 

4.5.9.  Final  Assignment  of  Photons  to  Channels 

The  final  assignment  of  each  photon  event  to  a  channel  is 
performed  by  cf_identify_channel,  just  as  in  §  4.5.2, 
but  with  two  modifications:  First,  we  consider  only  photon 
events  with  CHANNEL  >  0;  unassigned  events  (which  are  not 
motion  corrected)  remain  unassigned.  Second,  we  do  not  pad 
the  extraction  windows  by  an  additional  ±10  pixels  in  Y,  al¬ 
though  the  extraction  windows  for  histogram  data  are  padded 
by  ±  8  Y  pixels  (or  an  amount  equal  to  the  binning  factor  in 
Y,  if  other  than  8),  as  before. 

4.5.10.  Compute  Count  Rates 

For  time-tag  data,  cf_target_count_rate  computes 
the  count  rate  in  the  target  aperture  for  the  LiF  and  SiC  chan¬ 
nels.  To  account  for  dead-time  effects,  the  contents  of  the 
WEIGHT  array  are  used.  Events  in  airglow  regions  are  ex¬ 
cluded,  but  no  other  filters  are  applied  to  the  data.  Results  are 
written  to  the  LIF_CNT_RATE  and  SIC_CNT_RATE  arrays 
of  the  time-line  table.  For  histogram  data,  the  initial  values  of 
these  arrays,  taken  from  the  housekeeping  file  (§  A3),  are  scaled 
by  the  value  of  the  header  keyword  DET_DEAD. 

4.6.  Wavelength  Calibration 

Once  converted  to  the  FARF  and  corrected  for  optical  and 
spacecraft  drift,  the  data  can  be  wavelength  calibrated.  The 
module  cf_assign_wavelength  performs  three  tasks: 
first,  it  corrects  for  astigmatism  in  the  spectrograph  optics; 
second,  it  applies  a  wavelength  calibration  to  each  photon 
event;  and  third,  it  shifts  the  wavelengths  to  a  heliocentric 
reference  frame.  The  derivation  of  the  FUSE  wavelength  scale 
is  discussed  in  §  5.1. 

4.6.1.  Astigmatism  Correction 

The  astigmatic  height  of  FUSE  spectra  perpendicular  to  the 
dispersion  axis  is  significant  and  varies  as  a  function  of  wave¬ 
length  (Fig.  3).  Moreover,  spectral  features  show  considerable 
curvature,  especially  near  the  ends  of  the  detectors,  where  the 
astigmatism  is  greatest.  The  subroutine  cf_astigmatism 
shifts  each  photon  event  in  X  to  correct  for  the  spectral-line 
curvature  introduced  by  the  FUSE  optics,  providing  a  noticeable 
improvement  in  spectral  resolution  for  point  sources  (Fig.  6). 

The  astigmatism  correction  is  derived  from  observations  of 
GCRV  12336,  the  central  star  of  the  Dumbbell  Nebula  (M27), 
which  exhibits  H,  absorption  features  across  the  FUSE  wave 
band.  We  cross-correlate  and  combine  the  absorption  features 
from  a  small  range  in  X,  fit  a  parabola  to  each  set  of  combined 
features,  compute  the  shift  required  to  straighten  each  parabola. 


2007  PASR  119:527-555 


CalFUSE 


541 


Fig.  6. — Segment  of  LiF  1A  spectrum  before  ( bottom )  and  after  (top)  astigmatism  correction.  Note  the  reduction  of  curvature  in  the  absorption  features.  A 
detector  dead  spot  is  present  on  the  left  side  of  the  figure. 


and  interpolate  the  shifts  across  the  wave  band.  Because  an 
astigmatism  correction  has  been  derived  only  for  point  sources, 
no  correction  is  performed  on  the  spectra  of  extended  sources, 
airglow  spectra,  or  observations  with  APERTURE  =  RFPT. 

The  correction  is  stored  in  the  calibration  file  ASTG_CAL 
as  a  two-dimensional  image  representing  the  region  of  the  de¬ 
tector  containing  the  target  spectrum.  A  separate  image  is  pro¬ 
vided  for  each  aperture.  The  value  of  each  image  pixel  is  the 
astigmatism  correction  (AX  in  pixels)  to  be  applied  to  that  pixel. 
The  entire  image  is  shifted  in  Y  to  match  the  centroid  of  the 
target  spectrum,  and  the  appropriate  correction  is  applied  to 
each  photon  event  in  the  target  aperture.  The  corrected  X  co¬ 
ordinates  are  not  written  to  the  IDF,  but  are  passed  immediately 
to  the  wavelength-assignment  subroutine.  In  effect,  we  apply 
a  two-dimensional  wavelength  calibration  that  depends  on  both 
the  X  and  Y  coordinates  of  each  photon  event. 

4.6.2.  Assign  Wavelengths 

The  wavelength  calibration  is  stored  as  a  binary  table  ex¬ 
tension  in  the  calibration  file  WAVE_CAL  (§  5.1);  a  separate 
extension  is  provided  for  each  aperture.  Wavelengths  are  tab¬ 
ulated  for  integral  values  of  X,  assumed  to  be  in  motion-cor¬ 
rected  FARF  coordinates.  Given  the  astigmatism-corrected  X 
and  CHANNEL  arrays,  the  subroutine  cf_dispersion  con¬ 
siders  each  aperture  in  turn  and  reads  the  corresponding  cali¬ 
bration  table.  It  interpolates  between  tabulated  values  of  X  to 
derive  the  wavelength  of  each  photon  event.  Photon  events  not 
assigned  to  an  aperture  (CHANNEL  =  0)  are  not  wavelength 
calibrated. 

4.6.3.  Doppler  Correction 

The  component  of  the  spacecraft’s  orbital  velocity  in  the 
direction  of  the  target  is  stored  in  the  ORBITALJVEL  array 
of  the  time-line  table.  The  component  of  the  Earth’s  orbital 
velocity  in  the  direction  of  the  target  is  stored  in  the  IDF  file 
header  keyword  V_HELIO.  Their  sum  is  used  to  compute  a 
time-dependent  Doppler  correction,  which  is  applied  to  each 
photon  event  by  the  subroutine  cf_doppler_and_ 
heliocentric.  The  resulting  wavelength  scale  is  heliocen¬ 
tric.  Because  histogram  data  are  assigned  identical  arrival  times, 


their  Doppler  correction  is  not  time  dependent.  Histogram  ex¬ 
posures  are  kept  short  (approximately  500  s)  to  minimize  the 
resulting  loss  of  spectral  resolution.  The  final  wavelength  as¬ 
signed  to  each  photon  event  is  stored  in  the  LAMBDA  array 
of  the  IDF. 

4.7.  Flux  Calibration 

Because  the  instrument  sensitivity  varies  through  the  mis¬ 
sion,  we  employ  a  set  of  time-dependent  effective-area  files 
(AEFF_CAL).  We  interpolate  between  the  two  files  whose 
dates  bracket  the  exposure  start  time,  but  do  not  extrapolate 
beyond  the  most  recent  effective-area  file.  Within  each  cali¬ 
bration  hie,  the  instrumental  effective  area  is  stored  as  a  binary 
table  extension,  with  a  separate  extension  provided  for  each 
aperture.  The  module  cf_flux_calibrate  invokes  a  single 
subroutine,  cf_convert_to_ergs.  Considering  each  ap¬ 
erture  in  turn,  the  program  reads  the  appropriate  calibration 
hies,  interpolates  between  them  if  appropriate,  and  computes 
the  “flux  density”  of  each  photon  (in  units  of  ergs  cm-2),  ac¬ 
cording  to  the  formula 

ERGCM2  =  WEIGHT  x  /jc/LAMBDA/Aeff(X),  (1) 

where  ERGCM2,  WEIGHT,  and  LAMBDA  are  read  from  the 
photon-event  list  (§  A3),  h  is  Planck’s  constant,  c  the  speed 
of  light,  and  Aeff(A)  the  effective  area  at  the  wavelength  of 
interest.  Only  photon  events  assigned  to  an  aperture  are  hux 
calibrated;  events  with  CHANNEL  =  0  are  ignored.  The  hux 
density  computed  for  each  photon  event  is  stored  in  the 
ERGCM2  array  of  the  IDF.  This  array  is  not  used  by  the 
pipeline,  but  is  employed  by  some  of  our  interactive  IDF  ma¬ 
nipulation  tools. 

4.8.  Create  Bad-Pixel  Map 

When  possible,  spectra  are  extracted  using  an  optimal-ex¬ 
traction  algorithm  (§  4.9.4)  that  employs  a  bad-pixel  map 
(BPM)  to  correct  for  flux  lost  to  dead  spots  and  other  detector 
blemishes.  Because  motions  of  the  spacecraft  and  its  optical 
components  cause  FUSE  spectra  to  wander  on  the  detector,  a 
particular  spectral  feature  may  be  affected  by  a  dead  spot  for 
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Fig.  7. — Nighttime  scattered-light  image  for  detector  1A.  Note  the  vertical  scattered-light  stripe  to  the  right  of  the  image  center. 


only  a  fraction  of  an  exposure.  We  thus  generate  a  bad-pixel 
map  for  each  exposure.  The  module  cf_bad_pixels  reads 
a  list  of  dead  spots  from  a  calibration  file  (QUAL_CAL),  de¬ 
termines  which  of  them  overlap  the  target  aperture,  tracks  the 
motion  corrections  applied  to  each,  and  converts  the  motion- 
corrected  coordinates  to  wavelengths.  The  resulting  bad-pixel 
map  (identified  by  the  suffix  “bpm.fit”)  has  a  format  similar 
to  that  of  an  IDF  (§  A4),  but  the  WEIGHT  column,  whose 
values  range  from  0  to  1,  represents  the  fraction  of  the  exposure 
time  that  each  pixel  was  affected  by  a  dead  spot.  No  BPM  file 
is  created  if  the  (screened)  exposure  time  is  less  than  1  s.  BPM 
files  are  not  archived,  but  can  be  generated  from  an  IDF  and 
its  associated  jitter  file  using  software  distributed  with  CalFUSE 
(available  from  MAST). 

4.9.  Extract  Spectra 

Through  all  previous  steps  of  the  pipeline,  we  resist  the 
temptation  to  convert  the  photon-event  list  into  an  image.  In 
the  module  cf_extract_spectra,  we  relent.  Indeed,  we 
generate  four  sets  of  images:  a  background  model,  a  bad-pixel 
mask,  a  two-dimensional  probability  distribution  of  the  target 
flux,  and  a  spectral  image  for  each  extracted  spectrum  (LiF  and 
SiC).  Only  photon  events  that  pass  all  of  the  requested  screen¬ 
ing  steps  (§  4.4)  are  considered.  If  the  (screened)  exposure  time 
is  less  than  1  s  or  if  no  photon  events  survive  the  screening 
routines,  then  the  program  generates  a  null-valued  spectral  file. 

4.9.1.  Background  Model 

MicroChannel  plates  contribute  to  the  detector  background 
via  beta  decay  of  40K  in  the  MCP  glass.  On  orbit,  cosmic  rays 
add  to  this  intrinsic  background  to  yield  a  total  rate  of 
~0.5  counts  cm~2  s  '.  Scattered  light,  primarily  geocoronal 
Lya,  contributes  a  well-defined  illumination  pattern  (Fig.  7) 
that  varies  in  intensity  during  the  orbit,  with  detector-averaged 
count  rates  as  small  as  20%  of  the  intrinsic  background  during 
the  night,  and  1-3  times  the  intrinsic  rate  during  the  day.  We 
assume  that  the  observed  background  consists  of  three  inde¬ 
pendent  components:  a  spatially  uniform  dark  count  and  spa¬ 
tially  varying  day-  and  nighttime  scattered-light  components. 
Properly  scaling  them  to  the  data  is  thus  a  problem  with  three 


unknown  parameters.  We  attempt  to  fit  as  many  of  these  pa¬ 
rameters  as  possible  directly  from  the  data.  When  such  a  fit  is 
not  possible,  we  estimate  one  or  more  components  and  fit  the 
remainder. 

Background  events  due  to  the  detector  generally  have  pulse- 
height  values  lower  than  those  of  real  photons  (§  4.4.12).  The 
observed  dark  count  thus  depends  on  the  pulse-height  limits 
imposed  on  the  data.  An  initial  estimate  of  the  dark  count  as 
a  function  of  the  lower  pulse-height  threshold  is  read  from  the 
background  characterization  file  (BCHR_CAL).  The  day  and 
night  components  of  the  scattered-light  model  are  read  from 
separate  extensions  of  the  appropriate  time-dependent  back¬ 
ground  calibration  file  (BKGD_CAL).  The  background  models 
are  scaled  to  match  the  counts  observed  on  unilluminated 
regions  of  the  detector.  The  Y  limits  of  these  regions  (selected 
according  to  the  target  aperture)  are  read  from  header  key  words 
in  the  IDF.  Airglow  photons  in  these  regions  are  excluded  from 
the  analysis.  The  day  and  night  counts  in  the  background 
regions  of  the  detector  are  summed  and  recorded. 

In  its  default  mode,  the  subroutine  cf_scale_bkgd  es¬ 
timates  the  uniform  background  as  follows:  the  background 
regions  of  the  day  and  night  scattered-light  models  are  scaled 
by  their  relative  exposure  times,  summed,  and  projected  onto 
the  X  axis  (to  produce  a  histogram  in  X).  A  similar  histogram 
(called  the  “empirical  background  spectrum”)  is  constructed 
from  the  data.  An  iterative  process  is  used  to  determine  the 
uniform  component  and  scattered-light  scale  factor  that  best 
reproduce  the  observed  X  distribution  of  background  counts. 
The  uniform  component  is  then  subtracted  from  the  day  and 
night  totals  computed  above,  and  the  day  and  night  components 
of  the  scattered-light  model  are  scaled  to  match  the  remaining 
observed  counts. 

If  the  empirical  background  spectrum  is  too  faint,  we  do  not 
attempt  to  fit  it,  but  assume  that  the  uniform  component  of  the 
background  is  equal  to  the  tabulated  dark-count  rate  scaled  by 
the  exposure  time.  The  day  and  night  components  of  the  scat¬ 
tered-light  model  are  calculated  as  above.  Users  who  require 
a  more  accurate  background  model  may  wish  to  combine  data 
from  multiple  exposures  before  extracting  a  spectrum  (see 
§  6.6).  If  the  empirical  background  spectrum  is  very  bright — 
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as,  for  example,  when  nebular  emission  or  a  background  star 
contaminates  one  of  the  other  apertures  (and  thus  the  back- 
ground-sample  region) — no  fit  is  performed.  Instead,  the  day 
and  night  components  of  the  scattered-light  model  are  scaled 
by  the  day  and  night  exposure  times,  the  tabulated  dark-count 
rate  is  scaled  by  the  total  exposure  time,  and  the  three  com¬ 
ponents  are  summed  to  produce  a  background  image.  This 
scheme  is  also  used  for  histogram  data.  Because  histogram  files 
contain  data  only  from  the  region  about  the  desired  extraction 
window,  the  background  cannot  be  estimated  from  other 
regions  of  the  detector.  Fortunately,  histogram  observations 
typically  consist  of  short  exposures  of  bright  targets,  so  the 
background  is  comparatively  faint. 

Our  day  and  night  scattered-light  models  were  derived  from 
the  sum  of  many  deep  background  observations  spanning  hun¬ 
dreds  of  kiloseconds.  Individual  exposures  that  differ  markedly 
from  the  mean  were  excluded  from  the  sum.  The  data  were 
processed  only  through  the  FARF-conversion  and  data-screen- 
ing  steps  of  the  pipeline.  Airglow  features  were  replaced  with 
a  mean  background  interpolated  along  the  dispersion  (X)  axis 
of  the  detector.  An  estimate  of  the  uniform  background  com¬ 
ponent  was  subtracted  from  the  final  image,  and  the  data  were 
binned  by  16  detector  pixels  in  X.  This  process  was  performed 
on  both  day-  and  night-only  data  sets.  We  produce  a  new  set 
of  background  images  every  6  to  12  months,  as  the  effects  of 
gain  sag  and  adjustments  of  the  detector  high  voltage  slowly 
alter  the  relative  sensitivity  of  the  illuminated  and  background 
regions  of  the  detectors. 

Caveats :  While  early  versions  of  the  pipeline  (through  ver. 
2.4)  assume  a  10%  uncertainty  in  the  background  flux,  prop¬ 
agating  it  through  to  the  final  extracted  spectrum,  CalFUSE 
version  3  treats  uncertainties  in  the  background  as  systematic 
errors  and  does  not  include  them  in  the  (purely  statistical)  error 
bars  of  the  extracted  spectra. 

The  algorithm  assumes  that  the  intensity  of  the  uniform  back¬ 
ground  is  constant  throughout  an  exposure.  This  would  be  the 
case  if  it  were  due  only  to  the  detector  dark  count,  but  in  fact 
the  uniform  background  includes  a  substantial  contribution 
from  the  scattered  light  and  is  thus  brighter  during  daytime 
portions  of  an  exposure.  The  assumption  of  a  constant  uniform 
background  can  lead  the  algorithm  to  overestimate  the  scale 
factors  for  both  the  uniform  and  spatially  varying  components 
of  the  background  model.  A  better  scheme  would  be  to  fit  the 
day  and  night  components  of  the  uniform  background  sepa¬ 
rately.  Similarly,  when  the  empirical  background  spectrum  is 
very  faint,  adopting  the  tabulated  value  of  the  uniform  back¬ 
ground  is  not  the  best  solution.  It  is  likely  that  the  scattered- 
light  component  of  the  uniform  background  would  be  better 
estimated  from  the  observed  day  and  night  backgrounds.  The 
difference  between  the  tabulated  and  observed  levels  of  the 
uniform  background  will  become  greater  as  the  mission  extends 
into  solar  minimum  and  the  intensities  of  both  individual  air- 
glow  features  and  the  scattered-light  component  of  the  uniform 
background  continue  to  weaken. 


Grating  scattering  of  point-source  photons  along  the  disper¬ 
sion  direction  is  potentially  significant  and  is  not  corrected  by 
the  CalFUSE  pipeline.  Typical  values  are  1%-1.5%  of  the 
continuum  flux  in  the  SiC  channels,  and  10  times  less  in  the 
LiF  channels. 


4.9.2.  Probability  Array 

The  optimal-extraction  algorithm  (§  4.9.4)  requires  as  input 
a  two-dimensional  probability  array  representing  the  distribu¬ 
tion  of  flux  on  the  detector.  Separate  probability  arrays  derived 
from  high  signal-to-noise  ratio  stellar  observations  have  been 
computed  for  each  channel  and  stored  as  image  extensions  in 
the  weights  calibration  file  (WGTS_CAL). 

By  construction,  the  Y  dimension  of  the  probability  array 
represents  the  maximum  extent  of  the  extraction  window  for 
a  particular  aperture.  For  simplicity,  all  arrays  employed  by 
the  optimal-extraction  algorithm  are  trimmed  to  match  the 
probability  array  in  Y.  The  centroids  of  the  probability  dis¬ 
tribution  and  the  target  spectrum  (recorded  in  the  corre¬ 
sponding  file  headers)  are  used  to  determine  the  offset  be¬ 
tween  detector  and  probability  coordinates.  In  the  X 
dimension,  all  arrays  are  binned  to  the  output  wavelength 
scale  requested  by  the  user.  Default  wavelength  parameters 
for  each  aperture  are  specified  in  the  header  of  the  wave¬ 
length  calibration  file  (WAVE_CAL);  the  default  binning 
for  all  channels  is  0.013  A  per  output  spectral  bin,  which 
corresponds  to  approximately  2  detector  pixels.  The  back¬ 
ground  array,  originally  binned  by  16  pixels  in  X,  is  rescaled 
to  the  width  of  each  output  wavelength  bin  by  the  subroutine 
cf_rebin_background.  The  probability  array  is  re¬ 
scaled  by  the  subroutine  cf_rebin_probability_ 
array  to  have  a  sum  of  unity  in  the  Y  dimension  for  each 
wavelength  bin. 


4.9.3.  Bad-Pixel  Mask 

A  bad-pixel  mask  with  the  same  wavelength  scale  and  Y 
dimensions  as  the  probability  array  is  constructed  from  the 
BPM  file  (§  4.8)  by  the  subroutine  cf_make_mask.  The  array 
is  initialized  to  zero.  For  each  entry  in  the  BPM  file,  the  value 
of  the  WEIGHT  array  is  added  to  the  corresponding  pixel  of 
the  bad-pixel  mask.  The  mask  is  then  normalized  and  inverted 
so  that  the  center  of  the  deepest  dead  spot  has  a  value  of  0 
and  the  regions  outside  are  set  to  1 .  The  conversion  from  pixel 
to  wavelength  coordinates  can  open  gaps  in  the  mask,  which 
appear  as  values  of  unity  surrounded  by  pixels  with  lower 
values.  We  search  for  array  elements  that  are  larger  than  their 
neighbors  and  replace  them  with  the  mean  value  of  the  ad¬ 
joining  pixels.  If  the  BPM  file  is  absent  or  a  particular  aperture 
is  free  of  bad  pixels,  all  elements  of  the  bad-pixel  mask  are 
set  to  1. 
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4.9.4.  Optimal  (Weighted)  Spectral  Extraction 

The  extraction  subroutine,  cf_optimal_extraction, 
is  called  separately  for  each  of  the  two  target  spectra  (LiF  and 
SiC).  Inputs  include  the  photon-event  list  and  the  indices  of 
events  that  pass  through  the  various  screenings,  as  well  as  the 
two-dimensional  (2D)  background,  probability,  and  bad-pixel 
arrays  described  above.  For  numerical  simplicity,  extraction  is 
performed  using  the  WEIGHT  of  each  photon  event,  rather 
than  its  ERGCM2  value.  A  pair  of  2D  data  and  variance  arrays 
with  the  same  dimensions  as  the  probability  array  are  con¬ 
structed  from  the  good  photons  whose  CHANNEL  values  cor¬ 
respond  to  the  target  aperture.  For  time-tag  data,  this  process 
is  straightforward:  the  LAMBDA  and  Y  values  of  each  photon 
event  correspond  to  a  particular  cell  in  the  data  and  variance 
arrays.  That  cell  in  the  data  array  is  incremented  by  the  photon 
weight,  while  the  corresponding  cell  in  the  variance  array  is 
incremented  by  the  square  of  the  weight.  A  ID  raw-counts 
spectrum  (useful  for  the  statistical  analysis  of  low-count-rate 
data)  is  constructed  simultaneously:  for  each  photon  event 
added  to  the  data  array,  the  appropriate  bin  of  the  counts  spec¬ 
trum  is  incremented  by  1 . 

For  histogram  data,  the  process  is  more  complex,  because 
the  original  detector  image  is  generally  binned  by  8  pixels  in 
Y,  and  because  each  entry  in  the  photon-event  list  represents 
the  sum  of  many  individual  photons.  In  the  Y  dimension,  an 
event’s  WEIGHT  is  divided  among  8  pixels  (or  the  actual  Y 
binning  for  that  exposure,  if  different),  according  to  the  dis¬ 
tribution  predicted  by  the  probability  array.  In  the  X  dimension, 
each  event  is  assumed  to  have  a  width  in  wavelength  space 
equal  to  the  mean  dispersion  per  pixel  for  the  channel  (read 
from  the  DISPAPIX  keyword  of  the  WAVE_CAL  hie),  and 
the  WEIGHT  of  an  event  that  spans  the  boundary  between  two 
output  wavelength  bins  is  divided  between  them.  This  smooth¬ 
ing  in  X  helps  to  mitigate  the  “beating”  that  would  otherwise 
occur  between  detector  pixels  and  output  wavelength  bins. 

One-dimensional  background,  weights,  and  variance  spectra 
are  then  extracted  from  the  two-dimensional  background,  data, 
and  variance  arrays.  To  ensure  that  the  three  spectra  sample 
the  same  region  of  the  detector,  only  cells  in  the  2D  arrays  for 
which  the  corresponding  cell  in  the  probability  array  has  a  value 
greater  than  10  4  are  included  in  the  sum.  These  limits  differ 
slightly  from  those  defined  in  the  aperture  (CHID_CAL)  cal¬ 
ibration  files.  As  a  result,  the  ratio  of  the  final  weights  and 
counts  spectra  may  not  be  a  constant.  (Ideally,  their  ratio  would 
equal  the  mean  dead-time  correction  for  the  exposure.)  An 
initial  flux  spectrum  equal  to  the  difference  of  the  weights  and 
background  spectra  is  used  as  input  to  the  optimal-extraction 
algorithm. 

CalFUSE  employs  the  optimal-extraction  algorithm  de¬ 
scribed  by  Horne  (1986),  which  requires  as  input  the  2D  data, 
background,  probability,  and  bad-pixel  arrays  and  the  ID  initial 
flux  spectrum.  Originally  designed  for  CCD  spectroscopy,  the 
algorithm  has  been  modified  for  the  FUSE  detectors.  Specifi¬ 


cally,  instead  of  constructing  a  2D  spatial  profile  from  each 
data  set,  we  use  a  tabulated  probability  array;  the  2D  cosmic- 
ray  mask,  an  integer  array  in  the  original  algorithm,  is  replaced 
with  the  bad-pixel  mask,  which  is  a  floating-point  array;  and 
the  2D  variance  estimate  is  scaled  by  the  bad-pixel  mask.  Ex¬ 
traction  is  iterative:  in  the  original  version,  iteration  is  per¬ 
formed  until  the  cosmic-ray  mask  stops  changing.  In  our  ver¬ 
sion,  iteration  continues  until  the  output  flux  spectrum  changes 
by  less  than  0.01  counts  in  all  pixels.  If  the  loop  repeats  50 
times,  the  algorithm  fails.  The  number  of  iterations  performed 
is  written  to  the  OPT_EXTR  keyword  of  the  output-file  header. 

If  optimal  extraction  is  successful,  the  variance  of  the  optimal 
spectrum  is  computed  using  the  recipe  of  Horne  (1986).  We 
have  adapted  this  recipe  to  produce  weights  and  background 
spectra  such  that  FLUX  =  WEIGHTS  —  BKGD.  The  resulting 
background  spectrum  is  not  smooth.  Optimal  extraction  is  not 
performed  on  the  spectra  of  extended  sources  or  on  those  for 
which  the  quality  of  the  computed  spectral  centroid  is  not 
HIGH.  (Both  the  centroid  and  its  quality  flag  are  stored  in  file 
header  keywords.)  In  these  cases,  or  if  the  optimal-extraction 
algorithm  fails,  the  initial  flux,  variance,  weights,  and  back¬ 
ground  spectra  are  adopted. 

However  they  are  constructed,  the  final  FLUX,  ERROR 
(equal  to  the  square  root  of  the  variance),  WEIGHTS,  and 
BKGD  arrays  (all  in  units  of  counts)  are  returned  to  the  calling 
routine.  For  time-tag  data,  the  COUNTS  array  as  described 
above  is  returned.  For  histogram  data,  the  COUNTS  array  is 
computed  by  dividing  the  final  WEIGHTS  array  by  the  mean 
dead-time  correction,  which  is  stored  in  a  file  header  keyword 
TOT_DEAD.  Also  returned  is  the  QUALITY  array.  It  is  the 
product  of  the  probability  array  and  the  bad-pixel  map,  pro¬ 
jected  onto  the  wavelength  axis  and  expressed  as  an  integer 
between  0  and  100.  Its  value  is  0  if  all  the  flux  in  a  wavelength 
bin  is  lost  to  a  detector  dead  spot,  100  if  no  flux  is  lost. 

Caveats:  The  optimal-extraction  algorithm  is  designed  to 
improve  the  signal-to-noise  ratio  of  the  spectra  of  faint  point 
sources.  Unfortunately,  it  is  in  precisely  these  cases  that  the 
spectral  centroid  is  most  likely  to  be  uncertain.  Because  proper 
positioning  of  the  probability  array  is  essential  to  the  weighting 
scheme,  observers  of  faint  targets  may  wish  to  combine  the 
IDFs  from  multiple  exposures  and  recompute  their  spectral 
centroids  before  attempting  optimal  extraction. 

4.9.5.  Extracted  Spectral  Files 

Because  the  optimal-extraction  routine  returns  the  final 
FLUX  and  ERROR  arrays  in  units  of  counts,  the  spectral  ex¬ 
traction  module  applies  a  flux  calibration  to  both  arrays  using 
the  subroutine  cf_convert_to_ergs  described  in  §  4.7. 
Dividing  each  array  element  by  (EXPTIME  x  WPC),  where 
EXPTIME  is  the  length  in  seconds  of  the  (screened)  exposure 
and  WPC  is  the  width  in  angstroms  of  each  output  wavelength 
bin,  completes  the  conversion  to  units  of  ergs  cuT2  s  1  A~\ 
The  format  of  the  extracted  spectral  files  is  described  in  §  A5. 
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4.10.  Trailer  and  Image  Files 

A  number  of  supplementary  files  are  generated  by  CalFUSE 
and  archived  with  the  data.  For  each  exposure  and  detector 
segment,  the  pipeline  generates  a  trailer  file  and  a  pair  of  image 
files  in  Graphics  Interchange  Format  (GIF).  The  trailer  file 
(suffix  “.trl”)  contains  timing  information  for  all  pipeline  mod¬ 
ules,  in  addition  to  any  warning  or  error  messages  that  they 
may  have  generated.  The  first  image  file  contains  an  image  of 
the  detector  overlaid  by  a  wavelength  scale  and  extraction  win¬ 
dows  for  each  aperture  (suffix  “ext. gif’).  Only  photon  events 
flagged  as  good  are  included  in  the  plot,  unless  there  are  none, 
in  which  case  all  events  are  plotted.  The  second  image  file 
presents  count-rate  plots  for  both  the  LiF  and  SiC  target  ap¬ 
ertures  (suffix  “rat.gif ').  These  arrays  come  from  the  time-line 
table  in  the  IDF  and  exclude  photons  flagged  as  airglow.  These 
image  files  are  powerful  tools  for  diagnosing  problems  in  the 
data,  revealing,  for  example,  when  high  background  levels 
cause  the  SiC  1  LWRS  extraction  window  to  be  misplaced. 

4.11.  Observation-level  Files 

For  each  exposure,  CalFUSE  produces  LiF  and  SiC  spectra 
from  each  of  four  detector  segments,  for  a  total  of  eight  ex¬ 
tracted  spectral  files.  OPUS  combines  them  into  a  set  of  three 
observation-level  files  for  submission  to  MAST.  Observation- 
level  files  are  distinguished  from  exposure-level  files  by  having 
an  exposure  number  of  000.  Depending  on  the  target  and  the 
scientific  questions  at  hand,  these  files  may  be  of  sufficient 
fidelity  for  scientific  investigation.  Here  is  a  brief  description 
of  their  contents: 

ALL. — For  each  combination  of  detector  segment  and  chan¬ 
nel  (LiF  1  A,  SiC  1A,  etc.),  we  combine  data  from  all  exposures 
in  the  observation  into  a  single  spectrum.  If  the  individual 
spectra  are  bright  enough,  we  cross-correlate  and  shift  them 
before  combining.  (For  each  channel,  the  shift  calculated  for 
the  detector  segment  spanning  1000-1100  A  is  applied  to  the 
other  segment  as  well.)  If  the  spectra  are  too  faint  for  cross¬ 
correlation,  we  combine  the  individual  IDFs  and  extract  a  single 
spectrum  to  optimize  the  background  model.  Combined  spectra 
(WAVE,  FLUX,  and  ERROR  arrays)  for  each  of  the  eight 
channels  are  stored  in  separate  binary  table  extensions  in  the 
following  order:  1ALIF,  IB  LIF,  2BLIF,  2ALIF,  1ASIC,  1BSIC, 
2BSIC,  and  2ASIC. 

A  NO  (all,  night-only). — With  the  same  format  as  the  ALL 
files,  these  spectra  are  constructed  using  only  data  obtained 
during  the  nighttime  portion  of  each  exposure.  They  are  gen¬ 
erated  only  for  time-tag  data,  and  only  if  EXPNIGHT  >  0.  The 
shifts  calculated  for  the  ALL  files  are  applied  to  the  night-only 
data;  they  are  not  recomputed. 

NVO  (National  Virtual  Observatory). — These  files  contain  a 
single  spectrum  spanning  the  entire  FUSE  wavelength  range. 
The  spectrum  is  assembled  by  cutting  and  pasting  segments 
from  the  most  sensitive  channel  at  each  wavelength.  Segments 


are  shifted  to  match  the  guide  channel  (either  LiF  1  or  LiF  2) 
between  1045  and  1070  A.  Columns  are  WAVE,  FLUX,  and 
ERROR  and  are  stored  in  a  single  binary  table  extension. 

The  ALL  file  is  used  to  generate  a  “quick  look”  spectral  plot 
for  each  observation.  When  available,  combined  spectra  from 
channels  spanning  the  FUSE  wave  band  are  plotted  in  a  single 
GIF  image  file  (suffix  “specttagf.gif”  or  “spechistf.gif”).  This 
plot  appears  on  the  MAST  preview  page  of  each  observation. 
Four  additional  GIF  files  contain  the  combined  LiF  1,  LiF  2, 
SiC  1,  and  SiC  2  spectra  for  each  observation. 

Caveats :  Cross-correlation  may  fail,  even  for  the  spectra  of 
bright  stars,  if  they  lack  strong  spectral  features.  Examples  are 
nearby  white  dwarfs  with  weak  interstellar  absorption  lines.  If 
cross-correlation  fails  for  a  given  exposure,  that  exposure  is 
excluded  from  the  sum.  Thus,  the  exposure  time  for  a  particular 
segment  in  an  ALL  file  may  be  less  than  the  total  exposure 
time  for  that  observation. 

The  cataloging  software  used  by  MAST  requires  the  presence 
of  an  ALL  file  for  each  exposure,  not  just  for  the  entire  ob¬ 
servation.  We  generate  exposure-level  ALL  files,  but  they  con¬ 
tain  no  data,  only  a  FITS  file  header.  The  observation-level 
ALL  files  discussed  above  can  be  distinguished  by  the  string 
“OOOOOall”  in  their  names. 

4.12.  Quality  Control  and  Archiving 

Before  the  reduced  data  are  archived,  they  undergo  a  two- 
step  quality-control  process:  First,  a  set  of  automated  checks 
is  performed  on  each  exposure.  The  software  compares  the  flux 
observed  in  the  guide  channel  (LiF  1  or  LiF  2)  with  that  ex¬ 
pected  for  the  target  and  with  that  observed  in  the  other  three 
channels.  If  an  anomaly  is  detected,  a  flag  is  set  requesting 
manual  investigation.  The  software  works  well  for  bright  con¬ 
tinuum  sources,  but  often  flags  faint  or  emission-line  targets 
as  unsuccessful  observations.  Second,  a  member  of  the  FUSE 
operations  team  investigates  any  warnings  generated  by  the 
software.  If  it  is  determined  that  less  than  50%  of  the  requested 
data  were  obtained,  the  target  is  reobserved. 

The  philosophy  of  the  FUSE  project  is  to  archive  data  when¬ 
ever  possible,  even  if  they  do  not  satisfy  the  requirements  of 
the  original  investigator.  As  a  result,  the  MAST  archive  contains 
a  number  of  FUSE  data  sets  that  are  in  some  way  flawed  (e.g., 
misaligned  channels,  partial  loss  of  guiding,  or  no  good  ob¬ 
serving  time).  Users  should  be  aware  of  this  possibility. 

If  the  pipeline  detects  an  error  in  the  data  or  its  associated 
housekeeping  or  jitter  files,  it  writes  a  warning  to  both  the 
trailer  file  and  the  headers  of  the  IDF  and  extracted  spectral 
files.  Users  are  advised  to  scan  trailer  files  for  the  “WARNING” 
string,  and  spectral  files  for  “COMMENT”  records.  Occasion¬ 
ally,  the  FUSE  operations  team  inserts  comments  directly  into 
the  headers  of  raw  data  files.  Such  comments  may  warn  of  an 
unusual  instrument  configuration,  errors  in  the  reported  target 
coordinates,  or  data  obtained  during  slews. 

Observation  names  beginning  with  the  letter  “S”  are  science- 
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verification  observations  and  may  have  been  obtained  with  an 
unusual  instrument  configuration.  For  example,  the  program 
S523  was  designed  to  test  procedures  for  observing  a  bright 
object  by  defocusing  the  SiC  mirrors.  The  LiF  mirrors  were 
moved  for  some  S523  exposures.  As  a  result,  data  from  this 
program  should  be  used  with  caution.  Abstracts  for  all  FUSE 
observing  programs  are  available  from  MAST. 

5.  CALIBRATION  FILES 

5.1.  Wavelength  Calibration 

5.1.1.  Derivation  of  the  FUSE  Wavelength  Calibration 

Our  principal  wavelength  calibration  target  is  GCRV  12336, 
central  star  of  the  planetary  nebula  M27,  whose  spectrum  ex¬ 
hibits  a  myriad  of  molecular-hydrogen  absorption  features 
(McCandliss  et  al.  2007).  For  the  SiC  IB  channel,  these  data 
are  supplemented  at  the  shortest  wavelengths  by  spectra  of  the 
hot  white  dwarf  G  191-B2B  (Lemoine  et  al.  2002).  For  each 
channel,  an  M27  spectrum  was  fully  reduced,  corrected  for 
astigmatism,  and  used  to  derive  an  empirical  mapping  of  pixel 
to  wavelength.  A  theoretical  dispersion  solution,  derived  from 
standard  optical  expressions,  was  fitted  to  the  empirical  data 
with  only  its  constant  term  (the  zero  point  of  the  wavelength 
scale)  as  a  free  parameter.  This  shifted  theoretical  dispersion 
solution  was  then  used  to  generate  the  wavelength  calibration 
file. 

Early  versions  of  the  pipeline  relied  on  the  wavelength  cal¬ 
ibration  to  correct  for  nonlinearities  in  the  detector  X  scale  (the 
geometric  distortion  discussed  in  §  4.3.6).  Correcting  for  this 
effect  separately  has  greatly  improved  the  accuracy  of  the 
FUSE  wavelength  scale.  To  determine  the  geometric  distortion 
in  the  X  dimension,  a  spline  was  fitted  to  the  residuals  from 
each  aperture  (expected  minus  observed  X  coordinate  of  each 
absorption  feature;  Fig.  8).  In  practice,  residuals  from  all  six 
apertures  (both  LiF  and  SiC  channels)  were  included  in  the  fit, 
but  data  from  the  other  five  apertures  were  weighted  100  times 
less  than  those  for  the  aperture  being  fitted.  The  additional  data 
points  help  to  constrain  the  fit  in  wavelength  regions  where  the 
data  are  sparse  or  missing.  The  spline  fits  from  all  six  apertures 
were  then  used  to  construct  the  2D  map  of  detector  distortions 
in  the  X  dimension  that  is  used  by  the  geometric-distortion 
routine  described  in  §  4.3.6.  The  process  is  iterative,  with  the 
residuals  (ideally)  becoming  smaller  with  each  iteration. 

The  scatter  of  individual  measurements  about  the  spline  fit 
is  in  some  cases  caused  by  blended  absorption  lines,  and  in 
others  by  localized  distortions  induced  by  the  fiber-bundle 
structure  of  the  MCPs.  This  scatter  is  thus  a  fair  estimate  of 
the  inaccuracies  that  the  user  may  expect  in  the  relative  mea¬ 
surement  of  the  wavelength  of  any  given  feature.  The  wave¬ 
length  inaccuracies  caused  by  localized  distortions  are  3-4  de¬ 
tector  pixels  (0.025  A,  or  7  km  s  ')  at  most  wavelengths,  but 
may  be  as  large  as  6-8  pixels.  They  occur  in  tiny  windows 
about  1  to  3  A  wide,  depending  on  the  channel  and  segment. 
Some  data  sets  show  larger  residuals.  These  distortions  are 


Fig.  8. — Geometric  distortion  in  the  X-coordinate  of  the  LiF  1A  LWRS 
channel.  Data  represent  the  difference  between  the  measured  locations  of  H2 
lines  in  the  spectrum  of  GCRV  12336  and  those  predicted  by  a  theoretical 
dispersion  relation.  The  solid  line  is  a  spline  fit  to  the  residuals. 

inherent  in  the  FUSE  data  set  and  represent  the  ultimate  limit 
to  the  accuracy  of  the  FUSE  wavelength  calibration. 

5.1.2.  Zero-Point  Uncertainties 

The  FUSE  wavelength  calibration  assumes  that  the  motion- 
corrected  spectrum  falls  at  a  precise  location  on  the  detector. 
If  it  is  shifted  in  X,  then  the  wavelength  scale  of  the  extracted 
spectrum  will  suffer  a  zero-point  offset.  For  the  guide  channel 
(either  LiF  1  or  LiF  2),  the  dominant  source  of  wavelength 
errors  is  thermally  induced  rotations  of  the  spectrograph  grat¬ 
ings,  which  depend  on  the  satellite  attitude.  For  the  other  chan¬ 
nels,  additional  wavelength  errors  come  from  mirror  misalign¬ 
ments  that  shift  the  target  away  from  the  center  of  the  aperture. 
Such  misalignments  may  produce  zero-point  offsets  of  up  to 
±0.15  A  for  point  sources  in  the  LWRS  aperture.  Offsets  are 
less  than  ±  0.02  A  for  the  MDRS  aperture  and  are  negligible 
for  the  HIRS  aperture. 

We  define  the  zero  point  of  our  wavelength  scale  by  requiring 
that  the  Ly/3  airglow  feature,  observed  through  the  HIRS  ap¬ 
erture  and  processed  as  if  it  were  a  point  source,  be  at  rest  in 
spacecraft  coordinates  when  all  Doppler  corrections  are  turned 
off.  The  use  of  an  airglow  feature  eliminates  errors  due  to  mirror 
motions  in  the  nonguide  channels,  but  not  errors  in  the  grating- 
motion  correction,  so  we  measure  the  Ly/3  line  in  some  200 
background  exposures  and  shift  their  mean  velocity  to  0  km 
For  each  channel,  all  three  apertures  (HIRS,  MDRS,  and 
LWRS)  use  this  HIRS-derived  wavelength  scale  (WAVE_CAL 
ver.  022  and  greater). 

The  grating-motion  correction  (§  4.5.4)  is  designed  to  place 
the  centroid  of  each  Ly/3  airglow  feature  at  a  fixed  location  in 
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Fig.  9. — Measured  heliocentric  velocity  of  the  interstellar  O  I  X1039.23 
absorption  feature  in  each  of  47  exposures  of  the  white  dwarf  KPD  0005 +5106 
through  the  high-resolution  (HIRS)  aperture.  Holberg  et  al.  (1998)  report  a 
heliocentric  velocity  of  —7.50  ±  0.76  km  s-1  for  the  interstellar  features  along 
this  line  of  sight. 


FARF  coordinates.  On  average,  it  achieves  that  goal:  for  our 
sample  of  200  background  exposures,  the  measured  velocity 
of  the  Ly/3  line  has  a  standard  deviation  of  between  2  and 
3  km  s  ',  depending  on  the  channel.  Unfortunately,  some  com¬ 
binations  of  pole  and  beta  angle  are  not  well  corrected,  leading 
to  velocity  offsets  of  10  km  s  1  or  more,  and  additional  mo¬ 
tions — of  either  the  gratings  or  some  other  optical  component — 
can  shift  the  extracted  spectra  by  several  km  s  1  from  one 
exposure  to  the  next. 

Figure  9  presents  the  measured  wavelength  of  the  interstellar 
O  i  XI 039  absorption  feature  in  47  exposures  of  the  hot  white 
dwarf  KPD  0005  +  5106  obtained  through  the  HIRS  aperture. 
The  2001  and  2002  data  show  little  scatter  and  yield  a  mean 
velocity  of  —10.7  ±  1.9  km  s  '.  The  2003  data  (all  from  a 
single  observation)  span  nearly  10  km  s  '.  The  2004  data  (again 
from  a  single  observation)  are  tightly  correlated  but  offset  by 
~7  km  s  '.  These  data  were  obtained  at  a  spacecraft  orientation 
(beta  and  pole  angles)  that  is  generally  well  corrected  by  our 
grating-motion  algorithm;  apparently,  some  other  effect  is  at 
work.  The  2006  data  come  from  the  LiF  2  channel,  which 
became  the  default  guide  channel  in  2005  (§  6.1). 

We  do  not  recommend  the  general  use  of  airglow  lines  to 
fix  the  absolute  wavelength  scale  of  point-source  spectra,  for 
several  reasons:  First,  airglow  emission  fills  the  aperture,  so 
the  resulting  airglow  lines  provide  no  information  about  the 
position  of  the  target  relative  to  the  aperture  center.  Second, 
the  jitter  correction  (for  all  channels)  and  the  mirror-motion 
correction  (for  the  SiC  channels)  are  inappropriate  for  airglow 
emission.  Third,  the  Doppler  correction  for  the  spacecraft’s 
orbital  motion  can  degrade  their  resolution. 


TABLE  1 

Stellar  Parameters  Adopted  for  FUSE  Flux  Standards 


Tat  log  g  +,ad 

Name  (K)  (cm  s-2)  (km  s_1) 


GD  71  .  32,843  7.783  80.0 

GD  659  .  35,326  7.923  33.0 

GD  153  .  39,158  7.770  50.0 

HZ  43  .  50,515  7.964  20.6 

GD  246  .  53,000  7.865  -13.2 

G  191-B2B  .  61,200  7.5 


Note.— For  G  191-B2B,  we  use  the  model  employed  by  Kruk  et  al.  (1999) 
for  the  final  Astro-2  calibration  of  HUT. 


5.1.3.  Diffuse  Emission 

The  FUSE  wavelength  scale  is  derived  from  astigmatism- 
corrected  point-source  spectra.  Extended-source  (diffuse)  spec¬ 
tra  are  not  corrected  for  astigmatism.  If  point-source  data  are 
processed  with  the  astigmatism  correction  turned  off,  the  re¬ 
sulting  wavelength  errors  are  less  than  about  4  detector  pixels, 
consistent  with  the  uncertainties  in  the  wavelength  scale.  There¬ 
fore,  the  present  FUSE  wavelength  calibration  should  be  ad¬ 
equate  for  extended-source  spectra.  Airglow  lines  are  useful 
for  determining  the  wavelength  zero  point  for  extended  sources 
that  fill  the  aperture. 


5.2.  Flux  Calibration 

5.2.1.  Derivation  of  the  Effective-Area  Curve 

The  FUSE  flux  calibration  is  based  on  in-flight  observations 
of  the  well-studied  DA  (pure  hydrogen)  white  dwarfs  listed  in 
Table  1,  which  have  been  observed  at  regular  intervals  through¬ 
out  the  mission.  For  each  channel,  data  from  multiple  stars  are 
combined  to  track  changes  in  the  instrument  sensitivity,  using 
a  technique  similar  to  that  developed  by  Massa  &  Fitzpatrick 
(2000)  for  the  International  Ultraviolet  Explorer  ( IUE )  satel¬ 
lite.  The  algorithm  yields  a  series  of  time-  and  wavelength- 
dependent  sensitivity  curves,  as  well  as  the  spectrum  of  each 
star,  in  units  of  raw  counts,  as  it  would  have  appeared  on  a 
date  early  in  the  mission,  which  we  choose  to  be  T0  =  1999 
December  31.  (We  refer  to  the  latter  as  “T0  spectra.”) 

For  each  star,  we  generated  a  synthetic  spectrum  using  the 
programs  TLUSTY  (ver.  200)  and  SYNSPSEC  (ver.  48)  of 
Hubeny  &  Lanz  (1995).  The  non-LTE  pure-hydrogen  model 
atmospheres  were  computed  according  to  a  prescription  by 
I.  Hubeny  (2003,  private  communication),  using  200  atmo¬ 
spheric  layers  to  ensure  an  optimal  absolute  flux  accuracy.  The 
atmospheric  parameters  listed  in  Table  1 ,  consistent  with  HST, 
IUE.  and  optical  observations,  were  used  to  compute  the  mod¬ 
els  (J.  Holberg  2003,  private  communication).  For  G  191-B2B, 
we  used  the  model  employed  by  Kruk  et  al.  (1999)  for  the  final 
Astro-2  calibration  of  the  Hopkins  Ultraviolet  Telescope 
(HUT).  Observations  of  these  stars  with  the  Faint  Object  Spec¬ 
trograph  aboard  HST  have  shown  that  the  models,  including 
parameter  uncertainties,  are  consistent  to  within  2%  at  wave- 
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lengths  longer  than  Lycc  (Bohlin  et  al.  1995;  Bohlin  1996). 
Uncertainties  in  the  far-ultraviolet  wave  band  are  slightly 
higher,  as  discussed  by  Kruk  et  al.  (1999). 

For  each  channel,  the  effective  area  in  units  of  cm2  is  computed 
by  dividing  one  or  more  T0  spectra  in  units  of  counts  s  1  A-1 
by  a  synthetic  white  dwarf  spectrum  in  units  of  photons  cm~2 
s  1  A-1.  Sensitivity  curves  for  the  LiF  1A  and  SiC  1A  channels 
are  presented  in  Figure  10.  (Effective-area  curves  for  all  FUSE 
channels  are  available  from  MAST.)  The  sensitivity  of  the  LiF 
1A  channel  decreased  by  ~15%  over  the  first  3  years  of  the 
mission,  but  appears  to  have  stabilized;  that  of  the  SiC  1  channel 
has  declined  by  ~45%  since  launch  and  is  falling  still  (although 
slowly).  Effective-area  curves  (AEFF_CAL)  for  each  channel 
and  detector  segment  were  generated  at  3  month  intervals  until 
the  loss  of  the  third  reaction  wheel  in  2004  December;  we  plan 
to  generate  them  at  6  month  intervals  for  the  duration  of  the 
mission. 

Caveats:  We  do  not  attempt  to  correct  spectra  obtained 
through  the  MDRS  and  HIRS  apertures  for  changes  in  instru¬ 
ment  sensitivity,  but  employ  a  single  effective-area  curve  for 
each.  The  low  throughput  of  these  apertures,  combined  with 
the  likelihood  that  their  spectra  are  nonphotometric,  makes 
tracking  changes  in  their  sensitivity  both  more  difficult  and  less 
useful  than  for  the  LWRS  aperture. 

5.2.2.  Systematic  Uncertainties 

The  greatest  uncertainties  in  a  line  or  continuum  flux  derived 
from  a  FUSE  spectrum  are  due  to  systematic  effects.  An  es¬ 
timate  of  the  uncertainty  in  our  flux  calibration  can  be  obtained 
by  comparing  the  effective-area  curves  derived  from  different 
white  dwarf  stars.  Differences  among  the  curves  reflect  errors 
in  both  the  model  atmospheres  and  the  stellar  parameters  on 
which  they  are  based.  In  most  channels,  the  scatter  in  the  de¬ 
rived  effective  areas  is  between  2%  and  4%. 

The  photometric  accuracy  of  FUSE  spectra  is  subject  to 
numerous  effects  that  cannot  be  fully  corrected  by  the 
CalFUSE  pipeline.  A  target  centered  in  an  aperture  of  the 
guide  channel  (LiF  1  or  LiF  2)  may  not  be  centered  in  the 
corresponding  apertures  of  the  other  three  channels.  Since  the 
loss  of  the  first  two  reaction  wheels  in  2001,  spacecraft  drifts 
may  move  the  target  out  of  even  the  guide-channel  aperture. 
While  the  pipeline  does  attempt  to  flag  times  when  the  target 
is  out  of  the  aperture,  the  algorithm  used  is  conservative  in 
that  it  underestimates  the  time  lost  to  pointing  errors  (§  4.4.6). 
The  user  is  advised  to  consult  the  count-rate  plots  genera¬ 
ted  by  the  pipeline  (suffix  “rat. gif’;  §  4.10)  and  the 
LIF_CNT_RATE  and  SIC_CNT_RATE  arrays  of  the  IDF 
time-line  table  to  determine  the  photometric  quality  of  an  ex¬ 
posure.  Using  tools  available  from  MAST  or  the  user-defined 
good-time  intervals  discussed  in  §  4.4.7,  users  can  reject  time 
periods  when  the  count  rate  is  low,  or  rescale  the  flux  of  low- 
count-rate  exposures. 


Date  (Year) 


980  1000  1020  1040  1060  1080  1100 

Wavelength  (A) 


Fig.  10. — FUSE  sensitivity  as  a  function  of  time.  Top\  Effective  area  of  the 
LiF  1A  and  SiC  1A  channels,  averaged  over  the  wavelength  region  1030- 
1040  A.  The  gap  between  2004  October  and  2006  May  represents  the  period 
after  the  loss  of  the  third  reaction  wheel,  when  few  calibration  targets  were 
observed.  Bottom :  Effective-area  curves  for  the  LiF  1A  and  SiC  1 A  channels, 
dated  1999  and  2006.  (For  both  channels,  the  1999  curve  has  the  higher 
effective  area.) 

When  a  point-source  target  falls  near  the  top  or  bottom  edge 
of  an  aperture,  vignetting  in  the  spectrograph  may  attenuate 
the  target  flux  in  a  wavelength-dependent  way.  Astigmatism 
gives  FUSE  spectra  the  shape  of  a  bow  tie  (Fig.  3).  If  vignetting 
is  important,  then  the  spectrum  will  lie  below  the  center  of  the 
aperture  on  one  side  of  the  bow  tie,  and  above  it  on  the  other. 
Significant  flux  loss  is  possible  in  wavelength  regions  far  from 
the  center  of  the  bow  tie. 

Other  systematic  uncertainties  are  imposed  by  various  de¬ 
tector  flat-field  effects;  their  relative  importance  depends  on 
one’s  scientific  goals.  For  narrow  emission  lines,  flux  uncer¬ 
tainties  are  dominated  by  the  moire  pattern  (high-frequency 
ripples  due  to  beating  among  the  arrays  of  microchannel  pores 
in  the  MCP  stack;  §  6.4),  unless  the  observation  was  obtained 
using  an  FP  split  or  the  equivalent  was  achieved  via  grating 
and  mirror  motions.  For  broad  features,  the  moire  is  not  im- 
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Fig.  1 1 . — Point-source  spectra  showing  the  effects  of  the  worm.  Spectra  A 
and  B,  obtained  with  the  LiF  IB  channel,  show  deep  depressions  near  1145 
and  1 160  A,  respectively.  The  wavelength  of  maximum  attenuation  varies  with 
the  F-position  of  the  target  within  the  aperture.  Spectrum  C,  obtained  with  the 
LiF  2A  channel,  is  unattenuated. 

portant,  but  larger  scale  flat-field  features  are.  These  effects  are 
discussed  in  the  FUSE  Instrument  and  Data  Handbook.  Finally, 
when  fitting  a  spectral  energy  distribution,  the  greatest  uncer¬ 
tainty  is  caused  by  worms  (§  6.3),  which  may  depress  the 
observed  flux  over  tens  of  angstroms  by  50%  or  more. 

5.2.3.  Extended  Sources 

The  FUSE  flux  calibration  is  derived  from  point-source  tar¬ 
gets.  Because  the  distribution  of  flux  in  the  cross-dispersion 
direction  differs  for  point  and  extended  sources,  it  is  possible 
that  the  instrumental  sensitivity  may  also  differ;  this  question 
has  not  been  explored  in  detail.  Extended  spectra  are  less  af¬ 
fected  by  worms  (§  6.3)  than  are  point-source  spectra.  More¬ 
over,  because  the  spectrum  of  a  diffuse  emitter  is  spread  over 
a  larger  region  of  the  detector,  it  will  suffer  less  from  local 
flat-field  effects. 

6.  DISCUSSION 

6.1.  Spacecraft  Guiding  on  the  LiF  2  Channel 

The  switch  from  FES  A  to  FES  B  as  the  default  guide  camera 
in  2005  July  has  had  two  principal  effects  on  the  quality  of  FUSE 
data.  First,  tracking  with  FES  A  ensured  that  targets  remained  in 
the  center  of  the  LiF  1  aperture,  which  is  the  most  sensitive  channel 
in  the  astrophysically  important  1000-1 100  A  wave  band.  Track¬ 
ing  with  FES  B  will  keep  targets  centered  in  the  LiF  2  aperture, 
increasing  the  likelihood  of  data  loss  in  the  LiF  1  channel. 
Second,  in  order  to  optimize  the  optical  focus  of  FES  B,  the 
LiF  2  FPA  was  moved  out  of  the  focal  plane  of  the  LiF  2 
primary  mirror.  Observations  of  point  sources  with  the  LWRS 
aperture  are  unaffected,  and  the  point-source  spectral  resolution 
of  this  channel  is  unchanged,  but  the  throughput  of  the  narrow 
LiF  2  apertures  is  reduced.  The  effective  transmission  of  the 
apertures  has  not  been  characterized  in  detail,  but  is  approxi¬ 
mately  70%  for  LIF2  MDRS  and  15%  for  LiF  2  HIRS,  versus 
98%  and  60%  for  their  LiF  1  counterparts.  The  spectral  res¬ 


olution  for  diffuse  sources  is  expected  to  be  slightly  lower  in 
LiF  2  than  in  LiF  1. 

6.2.  Scattered  Solar  Emission 

In  addition  to  airglow  lines,  scattered  solar  emission  features 
are  present  in  the  SiC  channels  when  observing  at  high  beta  angles 
during  the  sunlit  portion  of  the  orbit.  Emission  from  C  hi  X977.0, 
Ly /3  X1025.7,  and  O  vi  XX1031.9,  1037.6  has  been  positively 
identified.  Emission  from  N  m  X991.6  and  N  ii  X1085.7  may 
also  be  present.  It  is  believed  that  sunlight  is  scattered  by  re¬ 
flective,  silver-coated  Teflon  blankets  lying  above  the  SiC  baf¬ 
fles.  At  low  beta  angles,  scattered  solar  emission  is  less  ap¬ 
parent,  because  the  blankets  are  shaded  by  the  SiC  baffles  and 
the  open  baffle  doors  and  because  the  radiation  strikes  the 
blankets  at  a  high  angle  of  incidence.  It  is  unknown  at  which 
beta  angle,  if  any,  the  solar  emission  completely  disappears. 
Because  the  LiF  channels  lie  on  the  shadowed  side  of  the 
spacecraft,  solar  emission  lines  are  not  seen  in  LiF  spectra. 
C  in  and  O  vi  emission  observed  in  the  SiC  channels  during 
orbital  day  should  always  be  compared  with  the  emission  ob¬ 
served  either  with  the  LiF  channel  or  during  the  nighttime 
portion  of  an  orbit. 

Since  the  failure  of  the  third  reaction  wheel  in  2004  Decem¬ 
ber,  FUSE  mission  controllers  have  experimented  with  the  use 
of  nonstandard  roll  angles  to  improve  spacecraft  stability.  These 
roll  angles  can  place  the  spacecraft  in  a  configuration  that 
greatly  increases  the  sunlight  scattered  into  one  of  the  SiC 
channels.  The  scattered  light,  mostly  Lyman  continuum  emis¬ 
sion,  appears  as  an  increase  in  the  background  at  wavelengths 
shorter  than  about  920  A;  strong,  resolved  Lyman  lines  are 
present  at  longer  wavelengths.  When  present,  it  is  generally 
seen  in  only  one  of  the  two  SiC  channels.  We  have  no  way  to 
model  or  subtract  this  emission. 

6.3.  The  Worm 

The  spectra  of  point-source  targets  occasionally  exhibit  a 
depression  in  flux  that  may  span  as  much  as  50  A  (Fig.  11). 
These  depressions  appear  in  detector  images  as  narrow  stripes 
roughly  parallel  to  the  dispersion  axis  (Fig.  12).  The  stripes, 
known  as  worms,  can  attenuate  as  much  as  50%  of  the  incident 
light  in  affected  portions  of  the  spectrum.  Worms  shift  in  the 
dispersion  direction  when  the  target  moves  in  the  aperture. 
They  are  due  to  an  unfortunate  interaction  between  the  hori¬ 
zontal  focus  of  the  spectrograph  and  the  innermost  wire  grid 
(the  quantum-efficiency  grid;  §  4.3.6).  Since  the  location  of 
this  focus  point  is  a  function  of  wavelength,  the  strength  of  a 
worm  is  exquisitely  sensitive  to  the  exact  position  of  the  spec¬ 
trum  on  the  detector.  We  cannot  determine  this  position  with 
sufficient  precision  to  reliably  correct  for  flux  lost  to  worms. 
Although  most  prominent  in  LiF  IB  LWRS  spectra,  worms 
can  appear  in  all  channels  and  apertures.  Observers  who  require 
absolute  spectrophotometry  should  carefully  examine  FUSE 
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Fig.  12. — Detector  images  showing  the  effects  of  the  worm.  In  these  negative  images,  worms  appear  as  bright  stripes  parallel  to  the  dispersion  axis.  The  data 
shown  correspond  to  spectra  A  and  B  in  Fig.  11  and  span  wavelengths  between  1134  and  1187  A. 


spectral  image  files  for  the  presence  of  worms.  The  redundant 
wavelength  coverage  of  the  various  FUSE  channels  can  be  used 
to  mitigate  their  effects. 

6.4.  The  Moire  Pattern  in  Histogram  Data 

Since  the  release  of  CalFUSE  version  3.0,  users  have  re¬ 
ported  strong,  non-Gaussian  noise  in  the  spectra  of  some  bright 
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Fig.  13. — Moire  pattern  in  the  LiF  2B  spectrum  of  the  star  HD  209339, 
obtained  in  histogram  mode.  The  associated  error  array  is  overplotted.  The 
moire  ripples  are  strongest  on  this  detector  segment,  but  are  also  seen  on 
segments  1A  and  IB. 


stars  observed  in  histogram  mode.  An  example  is  shown  in 
Figure  13.  The  high-frequency  ripples  have  a  period  of  ap¬ 
proximately  9  detector  pixels,  or  about  0.06  A.  These  ripples 
are  a  moire  pattern  due  to  beating  among  the  arrays  of  micro- 
channel  pores  in  the  three  layers  of  the  MCP  stack  (Tremsin 
et  al.  1999).  The  moire  fringes  are  strongest  on  segment  2B, 
but  are  also  visible  on  segments  1A  and  IB.  The  motion  cor¬ 
rections  applied  to  time-tag  data  tend  to  smooth  out  this  effect, 
but  it  can  be  quite  strong  in  histogram  data.  Where  it  is  present, 
users  are  advised  to  smooth  or  bin  their  spectra  by  at  least  one 
resolution  element  to  reduce  its  effects.  This  and  other  detector 
artifacts  are  described  in  the  FUSE  Instrument  and  Data 
Handbook. 


6.5.  A  Note  about  Time 

The  FUSE  spacecraft  uses  Coordinated  Universal  Time 
(UTC).  The  spacecraft  clock  is  updated  periodically  from  the 
ground  using  a  procedure  that  corrects  for  the  signal  transit 
time  from  the  ground  station  to  the  spacecraft.  The  ground 
station  time  comes  from  GPS  satellites.  The  Instrument  Data 
System  receives  a  1  Hz  signal  from  the  spacecraft  that  is  used 
to  align  the  IDS  clock  with  the  spacecraft  clock  to  an  accuracy 
of  ±5  ms.  In  time-tag  mode,  the  IDS  typically  inserts  a  time 
stamp  into  the  data  stream  once  per  second,  but  can  insert  time 
stamps  as  frequently  as  125  times  per  second.  Unfortunately, 
the  binary  format  of  the  time  stamp  rounds  the  time  value  to 
the  nearest  1/128  of  a  second.  The  two  periods  beat  against 
one  another,  causing  the  loss  of  three  time  stamps  each  second. 
Additional  timing  uncertainties  due  to  delays  in  the  detector 
electronics  have  not  been  measured,  but  are  assumed  to  be  on 
the  order  of  a  few  milliseconds.  For  most  time-tag  observations, 
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for  which  time  stamps  are  recorded  only  once  per  second,  these 
effects  can  safely  be  ignored. 

Raw  time-tag  files  are  constructed  by  assigning  the  value  of 
the  most  recent  time  stamp,  in  units  of  seconds  from  the  ex¬ 
posure  start  time,  to  each  subsequent  photon  event.  The  fre¬ 
quency  of  these  time  markers  determines  the  temporal  reso¬ 
lution  of  the  data.  Photon-arrival  times  are  not  modified  by  the 
pipeline:  values  are  UTC  as  assigned  by  the  IDS.  In  particular, 
photon-arrival  times  are  not  converted  to  a  heliocentric  scale. 

6.6.  Combining  Data  from  Multiple  Exposures 

For  each  FUSE  observation,  OPUS  combines  data  from  in¬ 
dividual  exposures  into  a  set  of  observation-level  spectra,  as 
described  in  §  4.11.  While  these  files  are  sufficient  for  many 
projects,  other  projects  may  benefit  from  specialized  data  pro¬ 
cessing.  Here  are  some  points  to  keep  in  mind  when  combining 
FUSE  data  from  multiple  exposures:  For  bright  targets,  the 
goal  is  to  maximize  spectral  resolution,  so  it  is  important  to 
precisely  align  the  spectra  from  individual  exposures  before 
combining  them.  The  wavelength  zero  points  of  segments  A 
and  B  are  consistent  across  each  of  the  FUSE  detectors  (§  5.1), 
so  shifts  measured  for  one  detector  segment  can  safely  be  ap¬ 
plied  to  the  other.  For  observations  made  before  2005  July,  the 
LiF  1  spectrum  is  likely  to  have  the  most  accurate  wavelength 
scale,  so  it  serves  as  the  standard  for  the  other  three  channels. 
For  later  observations,  the  LiF  2  spectra  are  likely  to  be  the 
most  accurate.  A  procedure  to  cross-correlate  and  shift  spectra 


by  hand  is  described  in  the  FUSE  Data  Analysis  Cookbook. 
When  cross-correlating  the  spectra  of  point-source  targets,  it 
is  important  to  exclude  regions  contaminated  by  airglow  fea¬ 
tures,  as  their  motions  are  unlikely  to  track  those  of  the  target. 
For  faint  targets,  the  goal  is  to  optimize  the  fidelity  of  the 
background  model  by  maximizing  the  signal-to-noise  ratio  on 
background  regions  of  the  detector,  a  goal  achieved  by  com¬ 
bining  the  IDFs  from  multiple  exposures  before  extracting  the 
spectra.  A  variety  of  C-  and  IDL18-based  tools  to  perform  these 
and  other  data-analysis  tasks  have  been  generated  by  the  FUSE 
project.  Software  and  documentation  are  available  from  MAST. 
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18  IDL  is  a  registered  trademark  of  ITT  Corporation  for  its  Interactive  Data 
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APPENDIX  A 
FILE  FORMATS 


All  FUSE  data  are  stored  as  FITS  hies  (Hanisch  et  al.  2001) 
containing  one  or  more  header  +  data  units  (HDUs).  The  first 
is  called  the  primary  HDU  (or  HDU  1);  it  consists  of  a  header 
and  an  optional  A'-dimensional  image  array.  The  primary  HDU 
may  be  followed  by  any  number  of  additional  HDUs,  called 
“extensions.”  Each  extension  has  its  own  header  and  data  unit. 
FUSE  employs  two  types  of  extensions:  image  extensions  (a 
2D  array  of  pixels)  and  binary  table  extensions  (rows  and 
columns  of  data  in  binary  representation).  CalFUSE  uses  the 
CFITSIO19  subroutine  library  (Pence  1999)  to  read  and  write 
FITS  hies. 

Al.  Raw  Time-Tag  and  Histogram  Files 

FUSE  raw  data  hies  are  generated  by  OPUS  using  both  data 
downlinked  by  the  telescope  and  information  from  the  FUSE 
Mission  Planning  Database  (§  4.1).  Information  regarding  the 
target,  exposure  times,  instrument  configuration,  and  engi- 


19  See  http://heasarc.gsfc.nasa.gov/docs/software/fitsio/fitsio.html. 


neering  parameters  is  stored  in  a  series  of  header  keywords  in 
the  primary  HDU.  All  header  keywords  are  described  in  the 
FUSE  Instrument  and  Data  Handbook.  In  raw  time-tag  hies 
(Table  2),  the  primary  HDU  consists  of  a  header  only,  with  no 
associated  image  array.  HDU  2  contains  the  photon-event  list, 
with  arrival  time  (in  seconds  from  the  exposure  start  time),  raw 
detector  coordinates,  and  pulse  height  for  each  event  in  turn. 
HDU  3  lists  good-time  intervals  (GTIs)  calculated  by  OPUS. 
Raw  time-tag  hie  names  end  with  the  suffix  “ttagfraw.ht.”  They 
can  be  as  large  as  10-20  Mbytes  for  the  brightest  targets. 

The  data  in  raw  histogram  hies  (suffix  “histfraw.ht”)  are 
stored  as  a  series  of  image  extensions  (Table  3).  The  primary 
HDU  contains  the  same  header  keywords  as  time-tag  hies, 
along  with  a  small  (8  x  64  pixel)  image  called  the  spectral 
image  allocation  (SIA)  table.  The  SIA  table  is  used  to  map 
regions  of  the  detector  to  onboard  memory.  Each  element  in 
the  SIA  table  corresponds  to  a  2048  x  16  pixel  region  on  a 
detector  segment.  If  the  element  is  set  to  1,  the  photons  from 
the  corresponding  region  are  saved;  if  0,  they  are  discarded. 
Additional  image  extensions  follow,  each  containing  the  binned 
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TABLE  2 

Format  of  Raw  Time-Tag  Files 


Array  Name  Format  Description 

Primary  Header  Data  Unit  (HDU  1) 

Header  only.  Keywords  contain  exposure-specific  information. 

HDU  2:  Photon-Event  List 


TIME .  FLOAT  Photon  arrival  time  (seconds) 

X  .  SHORT  Raw  X-position  (0-16,383) 

Y  .  SHORT  Raw  T-position  (0-1023) 

PHA  .  BYTE  Pulse  height  (0-31) 


HDU  3:  Good-Time  Intervals 

START  .  DOLTBLE  GTI  start  time  (seconds) 

STOP .  DOUBLE  GTI  stop  time  (seconds) 


Note.— Times  are  relative  to  the  exposure  start  time  stored  in  the  header 
keyword  EXPSTART. 

image  of  some  region  of  the  detector;  these  regions  may  over¬ 
lap.  In  general,  science  data  are  binned  by  8  pixels  in  Y  and 
are  unbinned  in  X;  binning  factors  for  each  exposure  are  stored 
in  the  header  keywords  SPECBINX  and  SPECBINY.  While 
the  format  given  in  Table  3  is  standard,  any  number  of  image 
extensions  may  be  present  in  a  histogram  file.  Raw  histogram 
data  files  are  1—1.5  Mbytes  in  size. 

A2.  Housekeeping  and  Jitter  Files 

For  each  exposure,  a  single  housekeeping  file  is  generated 
by  OPUS  from  engineering  data  supplied  by  the  spacecraft 
(§  4.1).  Housekeeping  files  (suffix  “hskpf.fit”)  contain  62  ar¬ 
rays,  including  spacecraft  pointing  information,  detector  volt¬ 
age  levels,  and  various  counter  values,  in  a  single  binary  table 
extension.  Arrays  are  tabulated  once  per  second,  although  most 
parameters  are  updated  only  once  every  16  s.  Only  a  few  of 
the  housekeeping  arrays  are  employed  by  the  pipeline.  The 
detector  high  voltage  and  LiF,  SiC,  FEC,  and  AIC  counter 
arrays  are  used  to  populate  the  corresponding  arrays  in  the  IDF 
time-line  table  (§  A3). 

From  pointing  information  in  the  housekeeping  file,  OPUS 
derives  a  jitter  file  (suffix  “jitrf.fit”)  consisting  of  a  single  binary 
table  extension  with  four  columns:  TIME,  DX,  DY,  and 
TRKFLG.  The  time  refers  to  the  elapsed  time  (in  seconds) 
from  the  start  of  the  exposure.  Since  the  engineering  data  com¬ 
monly  begin  up  to  a  minute  before  the  exposure,  the  first  few 
entries  of  this  array  are  negative.  DX  and  DY  are  the  offsets 
along  the  X  (dispersion)  and  Y  (cross-dispersion)  directions  in 
arcseconds.  These  offsets  are  defined  relative  to  the  com¬ 
manded  position  of  the  telescope  (presumably  the  target  co¬ 
ordinates).  Finally,  TRKFLG  is  the  tracking-quality  flag.  Its 
value  is  —  1  if  the  spacecraft  is  not  tracking  properly,  0  if 
tracking  information  is  unavailable.  Values  between  1  and  5 
represent  increasing  levels  of  fidelity  for  DX  and  DY. 

Additional  details  regarding  the  contents  and  format  of  the 
housekeeping  and  jitter  files  are  provided  in  the  FUSE  Instru¬ 
ment  and  Data  Handbook. 


TABLE  3 

Format  of  Raw  Histogram  Files 


Image  Size3 

HDU  Contents  (binned  pixels) 

ib  .  SIA  table'  8  x  64 

2  .  SiC  spectral  image  (12-20)  x  16,384 

3  .  LiF  spectral  image  (12-20)  x  16,384 

4  .  Left  stim  pulse  2  x  2048 

5  .  Right  stim  pulse  2  x  2048 


Note.— While  this  table  describes  the  format  of  a  typical  raw  histogram 
file,  any  number  of  HDUs  are  allowed. 

“  Quoted  image  sizes  assume  the  standard  histogram  binning:  by  8  pixels 
in  Y,  unbinned  in  X.  Actual  binning  factors  are  given  in  the  primary  file  header. 
b  Header  keywords  of  HDU  1  contain  exposure-specific  information. 
c  The  SIA  table  describes  which  regions  of  the  detector  are  included  in  the 
file. 

A3.  Intermediate  Data  File  (IDF) 

The  IDF  (suffix  “idf.fit”)  contains  three  FITS  binary  table 
extensions;  their  contents  are  listed  in  Table  4.  The  file’s  pri¬ 
mary  header  data  unit  (HDU  1 )  is  copied  directly  from  the  raw 
data  file,  (For  histogram  data,  the  SIA  table  is  discarded.)  Var¬ 
ious  keywords  are  populated  by  the  initialization  routine 
(§  4.2)  and  by  subsequent  pipeline  modules.  The  first  binary 
table  extension  (HDU  2)  contains  the  photon  events  themselves. 
For  time-tag  data,  the  TIME,  XRAW,  YRAW,  and  PHA  arrays 
are  copied  from  the  raw  data  file,  and  the  WEIGHT  array  is 
initialized  to  1.0.  For  histogram  data,  each  image  pixel  is 
mapped  back  to  its  coordinates  on  the  full  detector,  which  are 
recorded  in  the  XRAW  and  YRAW  arrays.  The  WEIGHT  array 
is  initialized  to  the  number  of  photon  events  in  the  pixel.  Zero- 
valued  pixels  are  ignored.  Histogram  data  are  not  “unbinned.” 
Each  entry  of  the  TIME  array  is  set  to  the  midpoint  of  the 
exposure,  and  each  entry  of  the  PHA  array  to  20.  (Both  arrays 
are  subsequently  modified.) 

The  first  pipeline  module  (§  4.3)  corrects  for  various  detector 
effects;  it  scales  the  WEIGHT  array  to  correct  for  detector  dead 
time  and  populates  the  XFARF  and  YFARF  arrays.  (The  flight 
alignment  reference  frame  represents  the  output  of  an  ideal 
detector.)  Each  photon  event  is  assigned  to  one  of  six  aperture- 
channel  combinations  or  to  the  background  (Table  5),  and  a 
corresponding  code  is  written  to  the  CHANNEL  array  (§  4.5). 
After  corrections  for  mirror,  grating,  and  spacecraft  motion,  the 
photon’s  final  coordinates  are  recorded  in  the  X  and  Y  arrays. 
Although  they  are  floating-point  arrays,  XFARF,  YFARF,  X, 
and  Y  are  written  to  the  IDF  as  arrays  of  8  bit  integers  using 
the  FITS  TZERO  and  TSCALE  keywords.  This  process  ef¬ 
fectively  rounds  each  element  of  XFARF  and  X  to  the  nearest 
0.25  of  a  detector  pixel,  and  each  element  of  YFARF  and  Y 
to  the  nearest  0.1  of  a  detector  pixel. 

The  screening  routines  (§  4.4)  use  information  from  the  time¬ 
line  table  (described  below)  to  identify  photons  that  violate 
pulse-height  limits,  limb-angle  constraints,  etc.  “Bad”  photons 
are  not  deleted  from  the  IDF,  but  are  merely  flagged.  Flags  are 
stored  as  single  bits  in  an  8  bit  byte.  We  use  two  sets  of  flags: 
TIMEFLGS  for  time-dependent  effects,  and  LOC_FLGS  for 
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TABLE  4 

Format  of  Intermediate  Data  Files 
Array  Name  Format  Description 

Primary  Header  Data  Unit  (HDU  1) 

Header  only.  Keywords  contain  exposure-specific  information. 

HDU  2:  Photon-Event  List 


TIME .  FLOAT  Photon  arrival  time  (seconds) 

XRAW  .  SHORT  Raw  X-coordinate  (0-16,383) 

YRAW  .  SHORT  Raw  F-coordinate  (0-1023) 

PHA  .  BYTE  Pulse  height  (0-31) 

WEIGHT .  FLOAT  Photons  per  binned  pixel  for  HIST  data;  initially  1.0  for  TTAG  data 

XFARF .  FLOAT  X-coordinate  in  geometrically  corrected  frame 

YFARF .  FLOAT  F-coordinate  in  geometrically  corrected  frame 

X  .  FLOAT  X-coordinate  after  motion  corrections 

Y  .  FLOAT  F-coordinate  after  motion  corrections 

CHANNEL  .  BYTE  Aperture  +  channel  ID  for  the  photon  (Table  5) 

TIMEFLGS  .  BYTE  Time  flags  (Table  6) 

LOC_FLGS  .  BYTE  Location  flags  (Table  6) 

LAMBDA .  FLOAT  Wavelength  of  photon  (A) 

ERGCM2  .  FLOAT  Energy  density  of  photon  (ergs  cnU2) 


HDU  3:  Good-Time  Intervals 

START  .  DOUBLE  GTI  start  time  (seconds) 

STOP .  DOUBLE  GTI  stop  time  (seconds) 


HDU  4:  Time-Line  Table 

TIME .  FLOAT  Seconds  from  exposure  start  time 

STATUS_FLAGS  .  BYTE  Status  flags 

TIME_SUNRISE  .  SHORT  Seconds  since  sunrise 

TIME_SUNSET  .  SHORT  Seconds  since  sunset 

LIMB_ANGLE  .  FLOAT  Limb  angle  (deg) 

LONGITUDE  .  FLOAT  Spacecraft  longitude  (deg) 

LATITUDE  .  FLOAT  Spacecraft  latitude  (deg) 

ORBITAL_VEL  .  FLOAT  Component  of  spacecraft  velocity  in  direction  of  target  (km  s-1) 

HIGHJVOLTAGE  .  SHORT  Detector  high  voltage  (unitless) 

LIF_CNT_RATE  .  SHORT  LiF  count  rate  (counts  s_1) 

SIC_CNT_RATE  .  SHORT  SiC  count  rate  (counts  s-1) 

FEC_CNT_RATE .  FLOAT  FEC  count  rate  (counts  s_1) 

AIC_CNT_RATE  .  FLOAT  AIC  count  rate  (counts  s-1) 

BKGD_CNT_RATE  .  SHORT  Background  count  rate  (counts  s-1) 

YCENT_LIF  .  FLOAT  F-centroid  of  LiF  target  spectrum  (pixels) 

YCENT_SIC  .  FLOAT  F-centroid  of  SiC  target  spectrum  (pixels) 


Notes.— Times  are  relative  to  the  exposure  start  time  stored  in  the  header  keyword  EXPSTART.  To  conserve 
memory,  floating-point  values  are  stored  as  shorts  (using  the  FITS  TZERO  and  TSCALE  keywords),  except  for 
TIME,  WEIGHT,  LAMBDA,  and  ERGCM2,  which  remain  floats. 


location-dependent  effects  (Table  6).  For  each  bit,  a  value  of 
0  indicates  that  the  photon  is  “good,”  except  for  the  day/night 
flag,  for  which  0  =  night  and  1  =  day.  It  is  possible  to  modify 
these  flags  without  rerunning  the  pipeline.  For  example,  one 
could  exclude  daytime  photons  or  include  data  taken  close  to 
the  Earth  limb. 

The  LAMBDA  array  contains  the  heliocentric  wavelength 
assigned  to  each  photon  (§  4.6),  and  the  ERGCM2  array  records 
its  “energy  density”  in  units  of  ergs  cm~2  (§  4.7).  To  convert 
an  extracted  spectrum  to  units  of  flux,  one  must  divide  by  the 
exposure  time  and  the  width  of  an  output  spectral  bin. 

The  second  extension  (HDU  3)  is  a  list  of  GTls.  The  initial 
values  are  copied  from  the  raw  data  file,  but  they  are  modified 
by  the  pipeline  once  the  various  screening  routines  have  been 


run.  By  convention,  the  START  value  of  each  GTI  corresponds 
to  the  arrival  time  of  the  first  photon  in  that  interval.  The  STOP 
value  is  1  s  later  than  the  arrival  time  of  the  last  photon  in  that 
interval.  The  length  of  the  GTI  is  thus  STOP  —  START. 

The  third  extension  (HDU  4)  is  called  the  time-line  table.  It 


TABLE  5 

Aperture  Codes  for  IDF 
CHANNEL  Array 


Aperture 

LiF 

SiC 

HIRS  . 

1 

5 

MDRS  . 

2 

6 

LWRS  . 

3 

7 

Not  in  an  aperture  . 

0 

0 
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TABLE  6 

Bit  Codes  for  IDF  Time  and  Location  Flags 


Bit  Value 


Time  Flags 

8  .  User-defined  bad-time  interval 

7  .  Jitter  (target  out  of  aperture) 

6  .  Not  in  an  OPUS-defined  GTI,  or  photon  arrival  time  unknown 

5  .  Burst 

4  .  High  voltage  reduced 

3  .  SAA 

2  .  Limb  angle 

1  .  Day/Night  flag  (N  =  0,  D  =  1) 


Location  Flags 

8  .  Not  used 

7  .  Fill  data  (histogram  mode  only) 

6  .  Photon  in  bad-pixel  region 

5  .  Photon  pulse  height  out  of  range 

4  .  Right  stim  pulse 

3  .  Left  stim  pulse 

2  .  Airglow  feature 

1  .  Not  in  detector  active  area 


Note.— Flags  are  listed  in  order  from  most  to  least  significant  bit. 

contains  status  flags  and  spacecraft  and  detector  parameters 
used  by  the  pipeline.  An  entry  in  the  time-line  table  is  created 
for  each  second  of  the  exposure.  For  time-tag  data,  the  first 
entry  corresponds  to  the  time  of  the  first  photon  event,  and  the 
final  entry  to  the  time  of  the  final  photon  event  plus  1  s.  (Should 
an  exposure’s  photon-arrival  times  purport  to  exceed  55  ks,  we 
create  time-line  entries  only  for  each  second  in  the  good-time 
intervals.)  For  histogram  data,  the  first  element  of  the  TIME 
array  is  set  to  zero,  and  the  final  element  to  EXPTIME+1 
(where  EXPTIME  is  the  exposure  duration  computed  by 
OPUS).  Because  we  require  that  EXPTIME  equal  both  2 
(STOP  —  START),  summed  over  all  entries  in  the  GTI  table, 
and  the  number  of  good  times  in  the  time-line  table,  we  must 
flag  the  final  second  of  each  GTI  as  bad.  No  photons  are  as¬ 
sociated  with  the  STOP  time  of  a  GTI. 

Only  the  day/night  and  OPUS  flags  of  the  STATUS_FLAGS 
array  are  populated  when  the  IDF  is  created;  the  other  flags 
are  set  by  the  various  screening  routines  (§  4.4).  The  elements 
of  the  TIME_SUNSET,  TIME_SUNRISE,  LIMB_ANGLE, 
LONGITUDE,  LATITUDE,  and  ORBITALJVEL  arrays  are 
computed  from  the  orbital  elements  in  the  FUSE.TLE  file.  The 
HIGH_VOLTAGE  array  is  populated  with  values  from  the 
housekeeping  file.  The  LIF_CNT_RATE  and  SIC_CNT_RATE 
arrays  are  initially  populated  with  values  derived  from  the  LiF 
and  SiC  counter  arrays  in  the  housekeeping  file.  For  time-tag 
data,  these  arrays  are  eventually  updated  with  the  actual  count 
rates  within  the  target  aperture,  excluding  regions  contaminated 
by  airglow.  The  FEC_CNT_RATE  and  AIC_CNT_RATE,  de¬ 
scribed  in  §  4.3.2,  are  also  derived  from  counter  arrays  in  the 
housekeeping  file.  For  time-tag  data,  the  BKGD_CNT_RATE 
array  is  populated  by  the  burst-rejection  routine  (§  4.4.5)  and 


represents  the  count  rate  in  predefined  background  regions  of 
the  detector,  excluding  airglow  features.  The  array  is  not  pop¬ 
ulated  for  histogram  data.  The  YCENTJLIF  and  YCENT_SIC 
arrays  trace  the  centroid  of  the  target  spectra  with  time  before 
motion  corrections  are  applied.  These  two  arrays  are  not  used 
by  the  pipeline. 

Raw  time-tag  files  (§  Al)  employ  the  standard  FITS  binary 
table  format,  listing  TIME,  X,  Y,  and  PHA  for  each  photon 
event  in  turn.  The  intermediate  data  files  have  a  slightly  dif¬ 
ferent  format,  listing  all  of  the  photon  arrival  times,  then  the 
X  coordinates,  then  the  Y  coordinates.  Formally,  the  table  has 
only  one  row,  and  each  element  of  the  table  is  an  array.  (To 
use  the  STSDAS  terminology,  IDFs  are  written  as  3D  tables.) 
The  MDRFITS  function  from  the  IDL  Astronomy  User’s  Li¬ 
brary20  (Landsman  1993)  can  read  both  file  formats;  some  older 
FITS  readers  cannot.  Note  that  because  HDUs  2  and  4  of  the 
IDFs  contain  floating-point  arrays  stored  as  shorts  (using  the 
TZERO  and  TSCALE  keywords),  calls  to  MRDFITS  must 
include  the  keyword  parameter  FSCALE. 

A4.  Bad-Pixel  Maps  (BPM  Files) 

The  BPM  files  (suffix  “bpm.fit”;  §  4.8)  consist  of  a  single 
binary  table  extension.  Its  format  is  similar  to  that  of  the  IDF, 
but  it  contains  only  five  columns:  X,  Y,  CHANNEL,  WEIGHT, 
and  LAMBDA.  The  WEIGHT  column,  whose  values  range 
from  0  to  1,  represents  the  fraction  of  the  exposure  time  that 
each  pixel  was  affected  by  a  dead  spot.  The  BPM  files  are  not 
archived,  but  can  be  generated  from  the  IDF  and  jitter  file  using 
pipeline  software  available  from  MAST. 

A5.  Extracted  Spectral  Files 

Extracted  spectra  (suffix  “fcal.fit”;  §  4.9)  are  stored  in  a 
single  binary  table  extension.  Its  contents  are  presented  in 
Table  7.  Note  that  the  spectra  are  binned  in  wavelength.  The 
bin  size  can  be  set  by  the  user,  but  the  default  is  0.013  A, 
which  corresponds  to  about  2  detector  pixels,  or  about  one- 
fourth  of  a  spectral  resolution  element.  The  WAVE  array  re¬ 
cords  the  central  wavelength  of  each  spectral  bin.  For  time-tag 
data,  the  COUNTS  array  represents  the  total  of  all  (raw)  photon 
events  assigned  to  the  target  aperture.  For  histogram  data,  the 
COUNTS  array  is  simply  the  WEIGHTS  array  divided  by  the 
mean  dead-time  correction  for  the  exposure.  If  optimal  extrac¬ 
tion  is  performed,  the  values  of  the  FLUX,  ERROR, 
WEIGHTS,  and  BKGD  arrays  are  determined  by  that  algo¬ 
rithm.  As  a  result,  the  ratio  of  WEIGHTS  to  COUNTS  is  con¬ 
stant  only  for  histogram  data.  The  QUALITY  array  records  the 
percentage  of  the  extraction  window  containing  valid  data.  It 
is  100  if  no  bad  pixels  fell  within  the  wavelength  bin,  and  0 
if  the  entire  bin  was  lost  to  bad  pixels. 


20  See  http://idlastro.gsfc.nasa.gov. 
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TABLE  7 

Format  of  Extracted  Spectral  Files 


Array  Name  Format  Description 

Primary  Header  Data  Unit  (HDU  1) 
Header  only.  Keywords  contain  exposure-specific  information. 


HDU  2:  Extracted  Spectrum 

WAVE  .  FLOAT  Wavelength  (A) 

FLUX  .  FLOAT  Flux  (ergs  cm"2  s_1  A-1) 

ERROR  .  FLOAT  Gaussian  error  (ergs  cm-2  s-1  A-1) 

COUNTS  .  INT  Raw  counts  in  extraction  window 

WEIGHTS  .  FLOAT  Raw  counts  corrected  for  dead  time 

BKGD  .  FLOAT  Estimated  background  in  extraction  window  (counts) 

QUALITY  .  SHORT  Percentage  of  window  used  for  extraction  (0-100) 
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